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ABSTRACT 



We present high S/N UV spectra for eight quasars at z ~ 3 ob- 
tained with VLT/FORS. The spectra enable us to analyze in detail the 
strongest emission features in the rest-frame range 1400-2000 A of each 
source (Cm]A1909, Silll]A1892, AI111AI86O, SillA1814, CivA1549 and 
SilvA1397). Previous work indicates that a component of these lines 
is emitted in a region with well-defined properties i.e., a high density 
and low ionization emitting region). Flux ratios AI111AI86O/ Sim] A1892, 
CivA1549/A1iiiA1860, SiivA1397/ Sim]A1892, CivA1549/SiivA1397 and 
SiiiA1814/ Sim] A1892 for this region permit us to strongly constrain elec- 
tron density, ionization parameter and metallicity through the use of di- 
agnostic maps built from CLOUDY simulations. Reliable estimates of 
the product density times ionization parameter allow us to derive the 
radius of the broad line region r blr from the definition of the ionization 
parameter. The tblr estimate and the assumption of virialized motions in 
the line emitting gas yields an estimate for black hole mass. We compare 
our results with estimates obtained from the r blr ~ luminosity correla- 
tion customarily employed to estimate black hole masses of high redshift 
quasars. 

Subject headings: galaxies: active — galaxies: high-redshift — quasars: 
general — quasars: emission lines 
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1. Introduction 

1.1. Interpreting Quasar Spectra 

Measuring relevant physical parameters from the observed broad-line spectra 
of quasars is still an open challenge. Identification and intensity measurements 
of the strongest emission lines has made possible a rough inference of the typical 
conditions of the emitting gas since the earliest days of quasar spectroscopy. The 
very first quasars of intermediate redshift discovered in the 1960s showed a fairly 
high ionization spectrum, with prominent lines of CivA1549, and HellA1640 in 
addition to strong Balmer lines of the low-redshift quasars. Photoionization by 
the central continuum source was considered the preeminent heating mechanism of 
the emitting gas. Significant Cin]A1909 emission suggested electron densities (n e ) 
in the range 10 9 - 10 10 cm" 3 . The observed intensity ratio Cm] A1909/CivA1549 
indicated ionization parameter (defined by Eq. [T] later in this paper) values of the 
order of 1CT 1 . This photoionization scenario was successful in explaining at least 
some quasar optical and UV spectra (see the review by Davidson and Netzer 1979 
for a synopsis) . More recent work emphasized the existence of several problems with 
this simple scenario. Low ionization lines (LILs), and especially Fell are too strong 
to be explained by a photoionized region of moderate density and column density 
(see for example Dumont and Mathez 1981; Joly 1987; Collin-Souffrin et al. 1988; 
Dumont and Collin-Souffrin 1990). These authors stressed that the LILs required a 
denser, low-temperature environment. 

We unfortunately lack a simple well-defined diagnostic measure of physical 
conditions in the broad line region. One strategy for estimating electron density 
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in astrophysical sources involves using two emission lines of the same ion, and 
with similar energies above the ground level, but with different radiative transition 
probabilities Aj&. In practice one often chooses two lines from the same term 
where one forbidden/semi-forbidden transition is associated with another that 
is semi-formidden/permitted in order to ensure very different values of A^i (for 
example, [Sim] A1882 and Sim] A1892). This technique is not straightforwardly 
applicable to the broad lines of quasars, precisely because the lines are broad, 
suitable candidates are too closely spaced in wavelength, and density is at least an 
order of magnitudes higher than the critical density for the forbidden transitions 
used in spectra of planetary nebulae and Hll regions. In addition the S/N and the 
resolution of quasar spectra are usually not very high. 

Feibelman and Aller (1987) used the Cm]A1909/ Sim] A1892 ratio to study 
the rather high-density environment typical of symbiotic stars (n e ~ 10 7 " 10 cm -3 ). 
In this case we have two semi-forbidden (intercombination) resonance lines with 
significantly different transition probabilities (see Table [T]). The lines are emitted 
by two iso-electronic species with somewhat different ionization potentials (11 eV 
for C + vs. 8 eV for Si + ). The ionic fraction is dependent on the relative abundance 
of silicon-to-carbon (to be assumed) as well as on the ionization structure within 
the emitting region (to be computed). Using lines of different ions introduces 
additional potentially serious sources of uncertainty. It is perhaps not surprising 
that most workers believe that n e in the BLR cannot be reliably estimated using 
quasar spectra. Strong Cin]A1909 emission would imply that n e cannot be very 
high (n e ~ 10 11_13 cm -3 ). Very high density was invoked to explain the rich low 
ionization spectrum (especially Fell) seen in the spectra of most quasars. Several 



-5 - 



lines in the UV spectrum of I Zw 1 point towards high density at least for the 
LIL emitting zone: prominent Fell, relatively strong AllIlA1860, and detection of 
ClllA1176 (Baldwin et al. 1996; Laor et al. 1997b). The region where these lines 
are produced cannot emit much Cin]A1909. But is Cin]A1909 really so strong in 
most quasars? BLR conditions are certainly complex and a single emitting region is 
not sufficient to explain both LILs and high ionization lines (HILs). 

1.2. Quasar Systematics 

Quasar spectra are not all alike. There are significant differences in line 
intensity ratios and broad line profiles from object to object (Bachev et al. 2004; 
Marziani et al. 2010). More importantly, these differences can be organized in a 
systematic way as has been realized since the early 1990s (Boroson and Green 1992). 
Since then several authors stressed the importance of the so-called eigenvector 1 
(El) of quasars (e.g. Gaskell et al. 1999). Sulentic et al. (2000, 2007) expanded 
the El trends into a 4-dimensional space involving optical, UV and X-ray measures. 
They also defined spectral types along a sequence occupied by AGN in an optical 
plane involving Fell and FWHM H/3 parameters. Objects at extreme ends of 
the El sequence are very different at almost all wavelengths and median spectra 
computed in spectral bins within this plane emphasize systematic changes in broad 
line properties (Sulentic et al. 2002, 2007). The differences have motivated the 
suggestion of a possible dichotomy between Narrow Line Seyfert 1 (NLSyls) like 
sources and broader line objects that include almost all radio-loud quasars. The 
most effective divider of the two quasar types appears to be at FWHM of the H/3 
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broad component H/?bc ~ 4000 km s 1 for low-to-moderate luminosity sources 
(Marziani et al. 2009). This corresponds to to Eddington ratio L/L^d ~ 0.2±0.1 
(Marziani et al. 2003b). 

The distinction between NLSyl-like objects (hereafter Population A sources 
with FWHM(H/3) < 4000 km s _1 ) and the rest of quasars (Population B with 
FWHM(H/3) > 4000 km s" 1 ) is of special relevance here. Pop. A sources show 
relatively low equivalent width lines with e.g. the A1900 blend ~30 A. A close 
analysis of the Cin]A1909 blend in the prototypical NLSyl I Zw 1 shows strong 
Sim] A1892, and AI111AI86O blended with rather weak Cni]A1909 along with 
prominent Felll(UV34)(A1895.5, A1914.0, and A1926.3) and Fell blends (Laor et al. 
1997b; Vestergaard and Wilkes 2001). The Clll]A1909 line is apparently so weak 
that FelllA1914 may be the most prominent feature at A ~ 1910A (cf. Hartig & 
Baldwin 1986). This interpretation is confirmed by detailed deblending of a source 
(SDSS J120144. 36+011611. 6) that can be considered a high luminosity analog of I 



Zw 1 (£7.1). 



Median composite UV spectra of low- z quasars show that AI111AI86O 
and Silll]A1892 are more prominent in 4DE1 spectral types A2 and A3. 
Bins Al, A2, A3 are defined in terms of increasing FeIIA4570 (see Fig. 1 
of Sulentic et al. 2002). I Zw 1, although belonging to the extreme type 
A3 (Bachev et al. 2004) is not unique as a NLSyl, since the fraction of 
Pop. A sources in bin A3 and A4 is « 20 % of all Pop. A sources in the 
sample of Zamfir et al. (2010). A significant number of extreme NLSyl 
are not classified as quasars in SDSS and must be collected from the 



-7- 



galaxy catalog (Hu et al. 2008). 

1.3. Emission Line Diagnostics and BLR Properties 

Emission lines and line ratios are used in diagnostic maps to estimate 
temperatures and electronic densities in galactic and extragalactic photoionized 
regions. Examples include HH-regions and galaxies with HH-region nuclear spectra 
where electron densities are less than n e K- 10 4 cm -3 . This method has been 
successfully applied to the narrow line region (NLR) in AGN. Application to the 
broad line quasars has yielded results that are difficult to interpret. One recent 
exception is the work by Maksuoka et al. (2008) who succeeded in analyzing the 
partly ionized region thought to emit most of the LILs in quasar spectra. They 
suggest that Ol 8446 and the Call triplet are emitted by dense, low ionization gas 
probably located in the periphery of the BLR. If the electron density and ionization 
conditions are known it is possible to derive, with additional assumptions, the 
distance of the BLR emitting region from the central continuum source (as stressed 
earlier also by Baldwin et al. 1996). 

The physical conditions of photoionized gas can be described by electron 
density n e , hydrogen column density N c , metallicity (normalized to solar), shape 
of the ionizing continuum, and the ionization parameter U. The latter represents 
the dimensionless ratio of the number of ionizing photons and the electron density 
n e or, equivalently, the total number density of hydrogen rin, ionized and neutral^] 

*In a fully ionized medium n w 1.2 71h- Wc prefer to adopt the definition based on nn because 
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Both U and are related through the equation 

[ + °° ¥dv 

where L v is the specific luminosity per unit frequency, h is the Planck constant, u 
the Rydberg frequency, c the speed of light, and r can be interpreted as the distance 
between the central source of ionizing radiation and the line emitting region. Note 
that Urine is the ionizing photon flux 

•<"> - ^ 0) 

If we know the product of tih and U, we can estimate the radius r of the BLR 
from Eq. [1} The dependence of U on tblr was used by Padovani Sz Rafanelli 
(1988) to derive central black hole masses assuming a plausible average value of the 
product tik-U. The typical value of n e was derived at that time from semiforbidden 
line Clll]A1909 which implied that the density could not be much higher than 
n e ~ 10 9 ' 5 cm -3 (Osterbrock & Ferland 2006). Padovani (1988) derived an average 
value < U- n e >m 10 9 ' 8 from several sources where tblr had been determined from 
reverberation mapping, and for which the number of ionizing photons could be 
measured from multiwavelength observations. The average value was then used to 
compute black hole masses for a much larger sample of Seyfert 1 galaxies and low- z 
quasars (Padovani & Rafanelli 1988; Padovani, Burg & Edelson 1989). Wandel, 
Peterson & Malkan (1999) compared the results of the photoionization method with 
the ones obtained through reverberation mapping, found a very good correlation 



it is the one employed in the CLOUDY computations. 
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for the masses computed with the two methods, and concluded that "both methods 
measure the mass of the central black hole." 



1.4. Outline of the paper 

The importance of the product U-nn goes beyond knowledge of the physical 
conditions in the BLR if it can lead to estimates of BLR radius and black hole 
mass. This paper identifies suitable emission line ratios that overcome some of 
the major problems in the analysis of emission lines of quasars. It also defines a 
photoionization method that can be applied to even the highest redshift quasars 
making use of high S/N near-IR spectroscopic data. In £j2] we present the spectra 
of 8 pilot sources obtained with the VLT/FORS; in §[3] we discuss data reduction; 
In §|4] we describe our method of fitting broad emission line profiles; in ^5] we give a 
phenomenological interpretation of the profile fits; in §[6] we describe a method for 
deriving BLR physical conditions and give the results of our fits; in §[7] we discuss 
two sources not belonging to our sample that show extreme behavior and that are 
helpful to understand more common quasar spectra; in £j8] we give the results for the 
photoionization method; in §[9] we derive the radius of the BLR (its distance from 
the ionizing source) and the mass of the black hole for each quasar in our sample; in 



S 10 we discuss our results; finally, in £ 11 we summarize our results and the prospect 



of a more extended application of our method. All the computations were made 
considering H =70 km s -1 Mpc -1 and a relative energy density = 0.70 and 
tt M = 0.3. 



2. Observations 



Data were obtained between Nov. 2006 and Jan. 2007 using the VLT2/FORS1 
telescope operated in service mode. FORS1 is the visual and near UV focal reducer 
and low dispersion spectrograph of the Very Large Telescope (VLT) operated by 
European Southern Observatory (ESO) (Appenzeller et al. 1998). Our VLT sample 
consists of 8 quasars with z ~ 3. In Figure [T] we show the spectra uncorrected for 
redshift. Tab. [2] provides a log of observations that is organized as follows. 

Column 1: object name, Col. 2: apparent B magnitude, Col. 3 redshift, Col. 4: 
line(s) used for redshift estimation: a) OlA1304.8, b) Cin]A1909; Col. 5: absolute B 
magnitude, Col. 6 flux at 6 cm taken from FIRST (Far InfraRed and Submillimetre 
Telescope), Col. 7: date (refers to time at start of exposure), Col. 8: Digital 
Integration Time, Col. 9: number of exposures with integration time equal to DIT, 
Col. 10: airmass at the beginning of each exposure, Col. 11 S/N in the continuum 
around 1700 A. 

The observation of one of our 8 quasars, J00521-1108, yielded only a low 
S/N spectrum which we retain because observed features in the blend at ~ 1900 
A are clear enough to fit the individual lines. Two sources, J01225+1339 and 
J02287+0002, are BAL quasars. We will keep them separate because CivA1549 is 
severely affected by absorption. 



3. Data Reduction 



Data were reduced using standard iraf tasks. All spectra were wavelength and 
flux calibrated in the observed frame and then corrected for Galactic extinction. 
Flux correction was applied using meteorological data provided by ESO. The 
observed flux was multiplied by the inverse of the light lost computed from the 
ratio seeing over slit width in arcsec. Correction to rest frame requires estimating 
the redshift z which is not a trivial task as outlined below. Rest frame correction 
also involved scaling the specific flux in flux per unit wavelength interval by a factor 
(1 + z) 3 . Measurements were carried out on the rest-frame spectra. It is necessity 
to describe below two important aspects of the data reduction. 



3.1. A &; B Atmospheric Bands Correction 

The A or B atmospheric band falls on top of the 1900A blend in many of the 
spectra. This is an important region for this study especially because it involves 
Sim] A1892, AI111AI86O, and SiiiA1814. In order to remove these absorption features 
we created an A+ B band template from standard star spectra used as specific flux 
calibrators. We scaled this template to find a best fit. Fig. [I] shows the A and B 
absorption correction where we make a line identification to illustrate which lines 
are affected. In cases where the A or B bands overlap a weak line like SillA1814 
the effect is considerable and measures of SiiiA1814 should not be considered at 
all or with extreme care. This happens for sources J00103-0037, J03036-0023, 
and J20497-0554. In cases where one of the bands overlaps a stronger line like 
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Sim] A1892 or AI111AI86O, the correction was good enough to permit accurate 
measures. 

3.2. Redshift Estimate 

Normally one uses strong narrow emission lines to set the rest frame in a 
quasar. In our case no strong narrow lines are available so we consider the peaks 
of Lya, CivA1549 and Clll]A1909. The Lya peak is affected by absorption and 
CivA1549 is a HIL feature often showing blueshifts and/or asymmetries (Gaskell 
1982, Espey et al. 1989, Corbin 1990, Tytler & fan 1992; Marziani et al. 1996; 
Richards et al. 2002, Baskin & Laor 2005; Sulentic et al. 2007). This is especially 
true in Pop. A sources. Clll]A1909 is blended with Sim] A1892 and Fein that is 
especially prominent in this region and could well affect the peak. Pop. B sources 
show a rather weak Fe spectrum making the Cm]A1909 peak a more reliable z 
estimator. 

Our best option is to use the low ionization line OlA1304 whenever it is strong. 
However it is blended with low ionization SillAA1304,1309 ( 2 Py 2 1/2 - 2 S 1/2 )- Both 
OlA1304 and SillAA1304,1309 are broad lines and in Pop. B sources might show 
large redshifts or even significant blueshifts. Simulations in the (n H , U) region of 
interest show OlA1304 2 SillA1304,1309 and this is confirmed in the spectrum I 
Zw 1 where OlA1304 and SillAA1304,1309 are resolved. The two components of the 
SillAA1304,1309 doublet are set to the same intensity (i.e., we assume an optically 
thick case). We model the blend OlA1304 + SillAA1304,1309 with 5 Gaussians; the 
three components of the Ol feature are produced by Bowen florescence mechanism, 
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and should show ratios consistent with their transition probabilities. Generating a 
model spectrum in iraf (lines broadened to 4000 km s" 1 ) yields a rest frame peak 
wavelength of 1304.8 ± 0.2 A (in vacuum) which we use as a reference for our VLT 
spectra assuming that there is no hint of systematic BC shifts as is the case for all 
Pop. A sources (the majority in our sample) and many Pop. B sources (Marziani 
et al. 2003a). 

Examination of Fig. [2]reveals that the peak of OlA1304 in source J00521-1108 is 
not observed clearly. We use Clll]A1909 to set the rest frame in this case. There are 
other sources J00103-0037, J02287+0002, J02390-0038 and J20497-0554 where the 
redshift estimation using both OlA1304 and Cm] A 1909 are not in good agreement. 
The largest disagreement was found for J02287+0002. Redshifts obtained for 
the three remaining quasars, J01225+1339, J03036-0023 and J23509-0052, were 
obtained from OlA1304. Fig. [3] shows the deredshifted VLT-FORS spectra for our 
sample of 8 quasars. 

4. Data Analysis 

4.1. Methodological Considerations on Multicomponent Fits 

The specfit IRAF task (Kriss 1994) allows us to fit the continuum, emission 
and absorption line components, Fell and Fein blends, etc. We fit two spectral 
ranges: (1) 1450-1680 A for analysis of CivA1549 and (2) 1750-2050 A for analysis 
of the 1900 A blend. Significant Fell and Fein emission are expected close to and 
underlying the 1900 A blend. Study of the 1900 A blend is especially difficult in 
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quasars because the lines are broad and the blending severe. We therefore need to 
take advantage of several previous results. 

4.2. Fell and Fen I Emission 

Our approach is completely empirical and employs an Fell + Fein template 
taken from templates successfully used in previous works. Our Fell template is 
based on a CLOUDY simulation and is not very far from the preferred model of 
Bruhweiler & Verner (2008). 

Fein lines are common and strong in the vicinity of Cin]A1909 as is evidenced 
by their presence in average LBQS (Francis et al. 1991) and SDSS (Vanden Berk et 
al. 2001) spectra. They appear to be strong when AllllA1860 is also strong (Hartig 
and Baldwin 1986). Vestergaard & Wilkes (2001) produced an Fein template based 
on the UV spectrum of IZwl. Since then, Sigut et al. (2004) have modeled the 
Fein BLR spectrum. See also Verner et al. (2003) for a plot of emission around 
1900 A. Lya pumping enhances Fein (UV 34)A1914.0 (Johansson et al. 2000) 
and this line can be a major contributor to the blend right on the red side of 
Clll]A1909 (see Fig. 2 of Vestergaard & Wilkes). We reproduced the option B of 
the empirical Fein template of Vestergaard & Wilkes (2001), taking advantage of 
the line identifications from Ekberg (1993). When detected we can use Fell UV 191 
to set a rough Fell level while the feature at 2080 A is helpful for a more precise 
estimation of the intensity of Fein. The continuum was fitted using the regions 
around 1450A (1750 and 1960A) that are relatively free of Fell emission (Vanden 
Berk et al. 2001). We used the same power-law to describe the continuum at both 
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the CIVA1549 and 1900A regions. 

4.3. Line Components 

We base our specfit analysis on several previous observational results. The 
most important ones are as follows: 

• Sulentic et al. (2002) gridded the broad component of FWHM B./3 B c versus 
i?Feii=W(FellA4570blend)/W(H/3 B c) parameter plane into bins of fixed A 
FWHM = 4000 km s -1 and A -Rf c ii= 0.5. Quasar spectra in different bins 
are different in many measures. As mentioned earlier, the largest differences 
are found between NLSyl-like objects, Pop. A, and broader sources of Pop. 
B with FWHM(H/3) > 4000 km s -1 . The gridding of Sulentic et al. (2002) is 
valid for low z (<0.7) quasars. At higher z an adjustment must be made since 
no sources with FWHM }i/3 BC < 3500 km s _1 exist above luminosity ~ 10 48 
ergs s -1 (Marziani et al. 2009). 

• Median spectra were computed for spectral bins from the atlas of Marziani 
et al. (2003a) who found that H/3 can be described by a Lorentz function in 
Pop. A sources and by the sum of 2 Gaussians in Pop. B sources (unshifted 
+ broader redshifted components) (Zamfir et al. 2010; Marziani et al. 2010). 

• A careful Fell subtraction reveals a blue-shifted H/3 component in some bin 
A3 sources (i.e. the stongest Fell emitters; Zamfir et al. 2010). 

• CivA1549 (HIL) and H/3 (L1L) profiles show significant differences in Pop. 
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A. Large CivA1549 blueshifts (< -1000 km s" 1 ) are observed in Pop. A 
only (Sulentic et al. 2007). HIL and LIL profiles are more similar in Pop. B 
sources. 

• We do not have H/3 observations for our high-z objects since there are no 
near IR spectra. We use the results of Marziani et al. (2003b, 2010): they 
show that the BC of Sim]A1892, AI111AI86O and CivA1549 lines is similar to 
the one of H/3, including the FWHM and profile shape, either Gaussian or 
Lorentzian. The similarity helps us to define whether an object is Population 
A or B in this paper. 

These observational results point toward three different components in broad 
line profiles (see Marziani et al. 2010) which can be described as follows: 

1. A broad component (BC) showing a roughly symmetric profile with FWHM 
in the range 1000-5000 km s -1 . It is consistent with the component identified 
by Matsuoka et al. (2008). This broad component dominates LILs in Balmer 
lines of Pop. A sources while it becomes less prominent in Pop. B. The profile 
is best modeled by a Lorentzian function in Pop. A sources while Pop. B 
profiles are better described by a Gaussian (Marziani et al. 2003b). 

2. A very broad component (VBC), as seen in LILs and HILs of most pop B 
sources but is absent from Pop. A profiles. The VBC can be modeled as a 
Gaussian (FWHM ~ 10000 km s _1 ) often with a significant shift to the red. 
It can be called a defining property of Pop. B sources (Marziani et al. 2010). 
This component is clearly identified in the CivA1549 line of Pop. B 
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objects, and is also appreciable on the red side of Clll]A1909 of Pop. 
B objects J00103-0037 and J02390-0038. 

3. A blueshifted broad component (BBC), denned as the residual emission in the 
CivA1549 line after subtracting a scaled BC Lorentzian profile (Marziani et 
al. 2010). This blueshifted component is often prominent in CivA1549 and 
Lya of Pop. A sources. It is much less intense in radio-loud Pop. B sources 
(Marziani et al. 1996; Punsly 2010; Richards et al. 2010). We model this 
profile as a blueshifted Gaussian. The Gaussian approximation is probably 
inappropriate especially if the BCC is strong: this component is believed to 
be produced in a partially-obscured radial flow, not in a virialized emitting 
system. In the present work we do not even try to fit the blueshifted 
in the doubly ionized lines. It is quite obvious from the fits that a 
possible contribution of this component would be negligible. 

Baldwin et al. (1996) presented a similar analysis. Their Fig. 2 
organizes spectra in a sequence that is roughly corresponding to El, 
going from AllllA1860-strong sources to objects whose spectra show 
prominent Cm]A1909 along with weak AlmA1860 (Bachev et al. 2004). 
Two of the three line components they isolated correspond to the ones 
we consider in this paper: a blue-shifted feature, and a more symmetric, 
unshifted and relatively narrow component that we call LIL-BC. Less 
obvious is the correspondence of a third feature, although it appears to 
be the redshifted part of what we call the VBC. 
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Several improvements have been introduced since the paper of 
Baldwin et al. 1996. These improvements are expected to make our 
analysis easier. First, the definition of a template of Fein emission 
(Vestergaard & Wilkes 2001), along with the possibility to model Fell in 
the 1400-2000 A spectral region with cloudy (Verner et al. 1999, 2004). 
The analysis of spectra along the El sequence allows one to see trends 
that make the interpretation of the emission line blends easier (Marziani 
et al. 2010). 

4.4. Expected emission from the various components 

We looked for evidence of three possible components as described 
above: BC, VBC and BBC only for the most intense HILs: CivA1549 
and SilvA1397. We did not include the contributions of the BBC for 
the doubly-ionized lines. In the case of Pop. A sources, we indeed 
assume that the BC contains the vast majority of the light. In the 
case of Pop. B sources we consider the contribution of the VBC. The 
AllllA1860 doublet shows no evidence of either a BBC or VBC. There 
is no obvious BBC for Silll]A1892 in the blend. There is evidence of a 
VBC of Clll]A1909 extending on the red side of the blend, and this is 
taken into account. Moreover, we expect that the AI111AI86O doublet is 
emitted exclusively in the BC, the region where Fell is also emitted. This 
is empirically confirmed by the aspect of the 1900 blend in many sources, 
where we do not see any evidence of BBC nor VBC. We remark that the 
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AllllA1860 doublet is relatively unblended, and that a BBC feature as 
strong as in CivA1549 would not easily escape visual detection. The same 
is also true for Silll]A1892. Several fits that included BBC in Clll]A1909 
yielded intensity, implying a large CivA1549/Ciii]A1909 (Marziani et 
al. 2010). The SilvA1397 + Oiv]A1402 blend closely resembles the shape 
of CivA1549, suggesting that BBC is relevant, especially for SilvA1397 
(Oiv]A1402 is expected to give a minority contribution to line emission at 
the high density derived for the BC; any Oiv]A1402 contribution to the 
BBC is not relevant to our method). BBC is very weak or undetectable 
in the vast majority of the H/3 profiles analyzed in Marziani et al. 2010 
(but see Zamfir et al. 2010 for several cases of H/3 BBC), while prominent 
in Lya; the Lya/H/3 ratio in this component is high. In summary, BBC 
is visually strong in Lya, CivA1549, and SilvA1397. HellA1640 in BBC 
is needed for a self-consistent fit of the ClvA1549+HellA1640 blend. The 
VBC of Silll]A1892 is poorly constrained, but in the fits where Clll]A1909 
VBC is visible, we always find Silll]A1892 < Clll]A1909, consistent with 
the high ionization level expected in the VBC region. Some SilvA1397 
VBC emission is assumed in the Pop. B fits. No VBC emission is 
expected in AI111AI86O and Fell. These constraints help also to make the 
fits less ambiguous. 

The absence of a VBC makes the decomposition of Pop. A spectra 
easier. As said earlier, the Fell intensity scale of the template (but Fell 
is in general weak) is anchored to the UV 191 intensity; similarly the 
Fein intensity is set by a feature external to the 1900 A blend (2080 
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A). There is no evidence of BBC of AI111AI86O and Silll]A1892; the 
AI111AI86O doublet profile is mostly unblended and defines the LIL-BC. 
The additional complication here is the Fein 1914 line whose intensity 
is affected by Lya pumping. Since the shift and FWHM are assumed 
the same for all lines (and templates) in the 1900 blend, the only free 
parameters in addition to shift and FWHM are the intensities of 6 
components, including the two from the templates, and SillA1814 and 
AI111AI86O that are not heavily blended. The specfit analysis is especially 
helpful to measure in a non-subjective way, taking all constraints into 
account, the two parameters that are most affected by the blend: the 
intensity of Silll]A1892 and Clll]A1909 (we repeat that any Fein A1914 
contribution in excess to the one of the adopted template is included in 
the estimated Clll]A1909 intensity). 

In the case of Pop B, the presence of a VBC does not really 
complicate the fit as a matter of fact. The Clll]A1909 line undeniably 
shows a VBC protruding on the red side of the 1900 blend. In any 
case, considering that we can expect the VBC to be assimilable to a 
shifted Gaussian with FWHM ~ 10000km s _1 , the unblended part of 
the Clll]A1909 VBC provides a strong constraint. The Silll]A1892 VBC 
is certainly the most difficult feature to ascertain, as it buried under 
Silll]A1892 and Clll]A1909 emission. We rely on the specfit results, that 
indicate negligible Silll]A1892 VBC. 



-21- 



4.5. Errors 

We identify five sources of error from the conditions for data reduction and 
methodological considerations described above: 



1. A & B atmospheric bands correction (already described in {3.1) 



2. Line profile shape, Gaussian or Lorentzian (Pop. A or B). The distinction 
between Pop. A and B is based on line width with the boundary at FWHM 
FL3 w km s _1 in low luminosity quasars and around ~ 5000 km s -1 at higher 
luminosity such as the eight sources presented here. Most of our quasars are 
unambiguously Pop. A or B because of line width and because Pop. B sources 
show an H/3 VBC and pop A sources a prominent CivA1549 BBC. In these 
cases only one profile shape (Gaussian or Lorentzian) was fitted. 

3. Rest-frame determination using OlA1304 or Cin]A1909. In some cases the 
redshift estimates derived from the two lines do not agree, most likely because 
of absorptions present in OlA1304.8 and because this is not a very intense line. 
The principal impact of uncertainty in the rest frame placement is estimation 
of the peak wavelength of CivA1549. If the line peak differs from A1549A, 
the BC intensity is diminished and we infer a greater contribution from 
the blue (BBC) component. Similarly, for the blend 1900A, the rest frame 
shift may increase or decrease our estimate for the strength of Clll]A1909 
with consequent decrease or increase of the Sim] A1892 contribution. This 
additional source of uncertainty affects J00103-0037, J02287+0002 and 
J20497-0554. 
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4. Fell intensity (continuum placement). Broad Fell emission can produce 
a pseudo-continuum affecting our estimates of emission line intensities. 
SiiiA1814 is especially affected in our spectra because it is weak. AI111AI86O is 
similarly affected when it is weak. The effect is less noticeable for CivA1549 
since expected Fell emission underlying the CivA1549 line is weak also for 
strong Fell emitters. 

5. Broad absorption lines (BALs) in quasars principally affect the blue side 
of CivA1549. We also find an absorption feature between FeIIA1787 and 
SillA1814 (eg. Fig. §). In sources J01225+1339 and J02287+0002 we can only 
set upper limits for line intensities (since we fit unabsorbed components). 

There are other sources of error such as small BC shifts and FWHM variations. 
We assume that the BC of the Silll]A1892, AI111AI86O, SillA1814 and 
CivA1549 line has the same FWHM and wavelength shift, although we 
allow for variations in their relative flux strengths and adopted profile 
type (Lorentzian for Pop. A and Gaussian for Pop. B). For Felll(UV34) 
and Fell we slightly relax this constraint. The Felll(UV34) and Fell 
emission is not very strong and the FWHM of individual features is 
poorly constrained by specfit. In the case of Clll]A1909 we need to 
consider the possibility that the profile is narrower because there might 
be a contribution from different regions: indeed, the SPECFIT routine 
usually converges toward a narrower profile if the Clll]A1909 width is 
not constrained. The effect depends on the strength of the Fein A 1914 
feature. 
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The goal of this study is to estimate diagnostic line ratios. A posteriori, we can 
say that the estimated line ratios are rather insensitive to the emission component 
profile shape: assuming a Gaussian or Lorentzian profile yields the same ratio for 
the strongest lines (i.e., CivA1549, Sim] A1892, AI111AI86O) upon which our analysis 
is based (with an uncertainty of ~10%). The same conclusion applies to the redshift 
uncertainty. Only in the case of J02287+0002 the redshift difference produces a 
significant effect due to a Az rs 0.16. However, we adopt the fit based on OlA1304: 
the alternative fit produces line ratios that are of not obvious interpretation and at 
variance with respect to the other sources. The most serious sources of error remain 
effects of A/B band overlap for SillA1814 and the presence of a BAL for CivA1549. 

5. Results of Line Component Analysis on Individual Objects 

In Figures from [I] to [6] we show our best fits for the VLT sample taking into 
account the considerations described in §[4| The fits of SilvA1397 are shown in Fig. 
[7] and the intensity values and equivalent widths are in Tables [3] and [5} We present 
here a phenomenological description of the fits. The line profiles, intensities 
and line ratios usually follow the trends of Pop. A or B sources, although 
in some objects there are features that are ambiguous (for example the 
line shape). In these cases however, we have assigned Pop. A or B type 
on the basis of the FWHM (see Cols. 5 and 6 of Table 
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5.1. Pop. A Objects 

• J03036-0023 - We estimate for this source a FWHM(BC)~3700 km s^ 1 and 
we use a Lorentzian function to fit the broad lines. The peak of CivA1549 
is blueshifted and requires a strong BBC (Fig. |4|a)). The bump on the red 
side of CivA1549 can be accounted for by HellA4686 BC and BBC. There is 
no evidence for a red shifted component in Cin]A1909 (Fig. |4|b)). AI111AI86O 
is prominent. Unfortunately the blue wing of SiiiA1814 and the red wing of 
FellA1787 are affected by A band absorption. 

• J20497-0554 - This source shows FWHM(BC)~3800 km s _1 . As for 
J03036-0023, the CivA1549 line can be accounted for by an unshifted BC 
(assumed Lorentzian) of a considerably contribution of a BBC (Fig. |4|c)). 
We see a prominent AI111AI86O line and FeIIA1787 (Fig. gjd)). Silll]A1892 
is affected by several narrow absorption lines; however, it is obviously strong. 
The lack of a red wing on Clll]A1909 suggests that no VBC is present. 

• J23509-0052 - This source has FWHM(BC)~3600 km s _1 . CivA1549 shows 
a slight blue asymmetry with a BBC (Fig. |4^e)) required to model it. The 
contribution of Fell is small and FellA1787 is weak (Fig. |4^f)). Clll]A1909 is 
very strong. AI111AI86O is affected by A band absorption; the profile we fit is 
probably an upper limit. This object could well belong to spectral type Al 
that includes Pop. A sources with the lowest Rpeiii ~ 0.5). 



5.2. Pop. B Objects 



J00103-0037 - This source has a FWHM(BC)~4500 kms" 1 . The red side 
of CivA1549 is blended with HellA4686. Fitting a BC with no shift plus 
a BBC to CivA1549 leaves a very large residual on the red side. A very 
broad redshifted component (VBC) is needed to model the spectrum (Fig|5] 
(a)). The faint narrow line under CivA1549 can be explained as the narrow 
component (NC) of ClvA1549(see Sulentic et al. 2007). The presence of 
a similar NC in Clll]A1909 could possible explain the large residual seen 
~1900. We specifically note the prominent Clll]A1909 emission and weak 
(but detected) AI111AI86O (Fig. gb)). The Fell "bump" at 1787 A(UV 
191) is appreciable. Fainter Fell emission is relatively unimportant because 
Fell creates a pseudo-continuum. SillA1814 is compromised by A-band 
absorption. The blend at 1900A includes a Clll]A1909 VBC and the fit 
indicates CivA1549/ Clll]A1909 (VBC) « 7 which is reasonable. 

J00521-1108 - This source shows the noisiest spectrum in the sample. We 
fit a FWHM(BC) ~ 5300 km s" 1 with CivA1549 requiring a large VBC to 
account for the red wing. (Fig. |5][c)). Absorption features seriously affect the 
ClvA1549 profile. The profile of Clll]A1909 is strongly asymmetric due to 
some sort of absorption on the red side. AI111AI86O is weak consistent with 
pop. B (Fig. gd)). 

J02390-0038 - This objects has a somewhat atypical Pop. B spectrum due to a 



-26 - 



very strong BBC in CivA1549 (Figure ge)). This source has a FWHM(BC)~ 
5400 km s -1 . Consistent with pop. B we find FellA1787 to be weaker than 
SillA1814. Cin]A1909 is flat topped and has a very similar intensity as 
Sim] A1892. Both CivA1549 and Clll]A1909 it show red wings indicating VBC 
emission. The A band absorption lies between Sim] A1892 and AI111AI86O 
(Fig. |5^f)) producing a possible overestimation of AI111AI86O. 

5.3. BAL QSOs 

• J01225+1339 - ClvA1549 is highly affected by two broad absorption lines 
(Fig. p[a)) with blueshifts of 5200 and 10800 km s _1 at peak absorption 
with equivalent widths/FWHM -12A / 3900 km s" 1 and -25A/5200 km s -1 , 
respectively. The blueshift of the CivA1549 peak leads us to suspect a large 
blueshifted BBC emission component. The 1900 A blend shows absorptions 
coincident with FellA1787 and the blue side of SillA1814 which is however 
unambiguously detected (Fig. plb)). AI111AI86O is prominent which implies 
that this a Pop. A source. The FWHM(BC) ~ 4400 km s" 1 is consistent 
with a high-luminosity Pop. A source. It is also possible that this BAL 
QSO is an outlier like Mark 231 at low-z (Sulentic et al. 2006), in 
other words an extreme Pop. A object. The Clll]A1909 is well fitted 
with a Lorentzian profile. Broad A band atmospheric absorption lies over 
Sim]A1892. 

• J02287+0002 - This object has a very complex spectrum. On one side it has 
a FWHM(BC) ss 4700 km s 1 which is typical of low-luminosity Pop. B 
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and the lines profiles are better fitted with Gaussians. On the other 
side, however, it shows features that are typical of extreme Pop. A 
sources: prominent FellA1787, strong AI111AI86O, no Clll]A1909 VBC 
(Fig. |6](d)). Considering that the FWHM limit between Pop. A 
and B is increasing with luminosity, the FWHM(BC) is within the 
limit of Pop. A. The Clll]A1909 line is not very flat topped but the 
similar intensities of Clll]A1909 and Silll]A1892 remind the case of 
J02390-0038. Also because it has a strong blue-shifted component 
in CIVA1549 atypical to Pop. B objects (Fig. [6(c)) . 

The estimated rest frame of this quasar differs by ~ 1400 km s -1 using 
OlA1304 and Clll]A1909. This is the largest discrepancy in our sample. In 
order to evaluate the effect of the z discrepancy we performed two fits using 
both rest frames. Figs. |6](a) and (b) use the Clll]A1909 restframe. In the 
1900 A blend, we found a contribution of Sim] A1892 similar to Clll]A1909. 
If we use the OlA1304 inferred rest frame we show in Figs. |6](e) and (f) that 
Clll]A1909 becomes stronger with a resultant decrease of Sim] A1892. A 
similar effect occurs for CivA1549 broad and blue-shifted components. 

The CivA1549 BAL shows a blueshift of 9100 km s _1 at deepest absorption, a 
EW of -14 A and a FWHM of 4600 km s _1 . 

Summing up, we are able to assign an A/B identification to all 
sources in our sample. The two BAL QSOs appear as objects of extreme 
Pop. A. We remind that the identification of Clll]A1909 in the BAL 
QSOs and in sources with strong AI111AI86O is debatable (Hartig and 
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Baldwin 1986): strong Fein 1914 could take the place of most Clll]A1909 
emission. 

6. Estimation of Physical Conditions in the Emitting regions 

6.1. CLOUDY Simulations 

We computed a multidimensional grid of CLOUDY (Ferland et al. 1998) 
simulations, (see also Korista et al. 1997) to derive U and nn from our spectral 
measurements. Simulations span the density range 7.00 < log %< 14.00, and 
—4.50 < XogU < 00.00, in intervals of 0.25. Each simulation was computed for 
a fixed ionization parameter and density assuming plane parallel geometry. The 
2D grid of simulations was repeated twice assuming N c = 10 23 and 10 24 cm -2 . 
Several cases were computed also for N c = 10 25 cm -2 . Metallicity was assumed to 
be either solar or five times solar. Two alternative input continua were used: 1) 
the standard AGN continuum of CLOUDY which is equivalent to the continuum 
described by Mathews and Ferland (1987) and 2) the low- z quasar continuum 
of Laor et al. (1997a). Computed line ratios are almost identical for fixed (U, 
nu). However the ionizing luminosity differs by more than a factor of 2 for a 
fixed specific continuum luminosity. The contour plots showing the distributions 
of Cm]A1909/ Sim]A1892, AI111AI86O/ Sim]A1892, SinA1814/ Sim]A1892, 
SiivA1397/ Sim]A1892, CivA1549/AlmA1860and CivA1549/ Sim]A1892 (Fig. 
[9j are generated from 29 x 19 = 551 simulations and assume a standard set of 
simulations using a Mathews and Ferland (1987) continuum, N c = 10 23 cm -2 , and 
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solar metallicity. The CLOUDY 08.00 computations included a model of the Fe + 



ion with 371 levels. The UV Fell template, described in §4.2| is based on a suitable 
CLOUDY simulation. Even if a relationship is very likely between the dense 
low ionization gas producing our diagnostic lines and Fell emission (supported 
observationally) our diagnostics do not use any Fell computation. The weak line 
OlA1304 is used only for rest frame estimation and not for diagnostic considerations. 

Apart from the hypothesis of plane-parallel geometry, no inferences are made 
about the actual distribution, location, and kinematics of the line emitting gas. 
The assumption of constant density is crude. If gas is distributed in clouds then 
magnetic confinement appears to be unnecessary to avoid cloud dispersion from 
pressure imbalance or cloud shear associated with a hot confining medium. Magnetic 
confinement could make density uniform within the cloud (Rees 1987; Bottorff and 
Ferland 2000). 



6.2. Intermediate Ionization Lines in the Blend at A1900 

The ratio Clll]A1909/ Sim] A1892 is density dependent because the transition 
probabilities of the two lines are so different: 114 s~ x vs 12600 s _1 (see Table Q. The 
forbidden lines at 1883 A and 1907A have such low transition probabilities that they 
are collisionally quenched at much lower density and will not be considered. Line 
ratios like Clll]A1909/ Sim] A1892 are useful diagnostics in a rather narrow range of 
density which depends on the transition probabilities. Above the critical density, 
emission lines originating from forbidden or semi-forbidden transitions become 
collisionally quenched, and hence weaker than lines for which collisional effects are 
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still negligible. Clll]A1909 is clearly unsuitable as a diagnostic for n e 3> 10 11 cm 3 , 
as the Cni]A1909/ Sim] A1892 0. Feldman et al. (1992) gives critical density 
of Sim] A1892 n e ~ 2 ■ 10 11 cm -3 . AI111AI86O is a permitted transition with large 
transition probability [A ~ 5 • 10 s s" 1 ) and has very-high critical density (i.e., its 
equivalent width goes to zero toward thermodynamic equilibrium, which occurs at 
very high density, when all line emission is collisionally quenched). Our 2D array of 
CLOUDY simulations shows that the ratio AI111AI86O/ Sim] A1892 is well suited 
to sample the density range 10 10 — 10 12 ' 5 cm -3 . Within this range the Sim] A1892 
intensity decreases smoothly by a factor 10; above the upper limit in density, the 
predicted intensity of Sim] A1892 decreases. This corresponds to the densest, low 
ionization emitting regions likely associated to the production of Fell. 

Intermediate-ionization lines such as Clll]A1909, AI111AI86O, Sim] A1892 and 
SilvA1397avoid the issue of collisional ionization (invoked for the production of 
low-ionization species by S. Collin and collaborators, as mentioned in the 
introduction). 

The widely separated doublet AI111AI86O is expected to be produced entirely 
in the fully ionized zone (Laor et al. 1997b). CLOUDY simulations confirm this 



suggestion. Figure 10 (top panel) shows the ionic fraction as a function of the 
geometrical depth in a cloud (or slab) of fixed column density (N c = 10 23 cm" 2 ) 
and density (hh= 10 12,5 cm -2 ). As it can be seen, Al ++ , Fe ++ , and Si ++ share a 
region of dominance deep in the cloud, close to the end of the Stromgren sphere. 
Beyond, in the partially-ionized zone (PIZ) there is, as by definition, a significant 
fraction of ionized hydrogen. The dominant ionic stages of Si and Fe become Fe + 
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and Si + . It is very appropriate to consider the AI111AI86O/ Sim] A1892 ratio since 
the two lines are apparently emitted in the very same zones within the gas slab. 
Available reverberation mapping results may support this interpretation but are 
rather difficult to extrapolate since they are limited to an handful of low luminosity 
objects. The main finding is that Clll]A1909 responds to continuum 
changes on timescales much longer than CivA1549 and other HILs. 
This results comes from the analysis of total Clll]A1909 + Silll]A1892 
in NGC 3783 (Onken & Peterson 2002), from the Clll]A1909 of NGC 
4151 (Metzroth et al. 2006) and of NGC 5548. It is intriguing that the 
Clll]A1909 cross-correlation delay in NGC 5548 (by far the best studied 
object) is even larger than that of H/3 (32 Id vs. 20 Id; Peterson et al. 
2002; Clavel et al. 1991). For fixed density, lines of higher ionization 
form at higher photon flux. The C + , Al + , Si + ionization potential are 
24, 18, and 16 eV respectively. These comparable ionization potentials 
are well below the one of HILs, X % _1 > 50 eV. However, the much lower 
Clll]A1909 critical density implies that the Clll]A1909 line should be 
formed farther out than Silll]A1892 and AI111AI86O if all these lines are 
produced under similar ionization conditions. 

6.3. The Sill Contribution to the UV Spectrum of Quasars 

AGN have a rather rich singly-ionized silicon spectrum in the range 1000-2000 
A, due to several resonant transitions from the ground term 3s 2 3p 12 P° to terms 
associated to the 3s l 3p 2 (A1814 and A1304), 3s 2 4s (A1531), 3s 2 3d (A1263) electronic 
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configurations. A feature on the red side of Lya, at 1263 A is detected in several 
Pop. A sources like I Zw 1 (Marziani et al. 2010). 

The A1309 feature is partially blended with the Ol triplet associated with 
a Bowen fluorescence mechanism from HI Ly/3. The A 1309 line is well resolved 
in sources like I Zw 1 with the extension of the broad A1309 feature suggesting 
significant Sill emission around A1309. The A1531 feature is blended with CivA1549. 
High S/N spectra of low redshift sources from HST (Laor et al. 1994) as well as for 
quasars at z ~ 4 (Constantin et al. 2002) indicate that the feature is very weak in 
most sources. 

The Sill feature at A1814 is detected in at least four of the five quasars 
studied by Laor et al. (1994). In the 6 objects used for the El sequence of 
Marziani et al. (2010), it is detected without reasonable doubt in I Zw 1 only. 
However this has to do more with the poor S/N of the 1900 A region than of 
anything else; good S/N HST spectra show an unambiguous detection in Akn 
120 and Mark 509, for example. Early detection of the UV Sill lines from IUE 
observations suggested collisional excitation and no relevant fluorescence effects 
(as observed in type la supernova?; Dumont and Mathez 1981). Fluorescence 
effects and recombination are revealed by optical lines which are emitted through 
several cascade routes to the ground state. However, a median spectrum of 
A2 quasars with S/N « 200 (Zamfir et al. 2010) barely detects any expected 
emission feature in the range 4000-6500 A. The prominence of the A1814 feature 
may then be associated to the relatively low temperature believed to exist in 
the innermost part of the line emitting "cloud" , T e ~ 5000 — 8000° K. We may 
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consider the equivalent width (W) ratio between the doublet lines at 1264 and 
at 1814 A both due to a 2 Ds — > 2 Pa transition: in the simplest case we have: 

2 2 

W 1263 /W 18U ~ (A u , 1264 /A u , 1814 ) (/ 1264 //i8i6) (Ai 26 4/Ai8i 6 ) 2 ~ 0.485^e"w « 0.27 
if T = 5000°, where / is the oscillator strength, N u are the density of the upper 
electronic levels giving rise to the two lines, and the energy difference between the 
transitions is AE 3 eV. 

6.4. Use of SilvA1397, CivA1549 and SillA1814 

The same caveats mentioned for the ratio Clll]A1909/ Sim] A1892 apply to the 
AllIlA1860/ Sim] A1892 ratio as well. The ground term 1 S has energy 16 eV and 
24 eV for Si +2 and C +2 respectively. A dependence on the ionization parameter 
is expected, as already mentioned. However, given the similarity in the ionization 
structure of the photo-ionized slab, it is after all not surprising that the ratio 
AI111AI86O/ Sim] A1892 is almost insensitive to the ionization parameter over a wide 
range of density. We remind that the detection of strong AI111AI86O alone already 
suggests that we are considering very high density emitting gas even if metallicity 
is super-solar. Simulations indicate that Aim A 1860 should increase smoothly with 
density and be weakest in canonical BLR if the density is 10 9 cm -3 (cf. Korista 
et al. 1997). 

The ratio AI111AI86O/ Sim] A1892 is therefore diagnosing high density gas, while 
the Cm] A 1909/ Sim] A1892 ratio covers the domain of n H ~ 10 10 cm" 3 . The ratio 
AI111AI86O/ Sim] A1892 alone is, generally speaking, insufficient to determine %. A 
second diagnostic ratio is needed to constrain U and to unambiguously derive nu- 
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We consider SillA1814/ Sim]A1892, CivA1549/ Silll]A1892, SilvA1397/ Silll]A1892 
as three diagnostic ratios suitable for constraining U. 

The SiiiA1814 doublet is conveniently placed, although weak or undetectable 
in most sources. The SiiiA1814/ Sim] A1892 ratio is anti-correlated with U in a 
regular form, as our CLOUDY simulations show, especially for log U < —2, and for 
density not much above the critical density of Sim] A1892. The behavior of the ratio 
SillA1814/ Sim] A1892 resembles the distribution of Cm] A1909/CivA1549 in the 
plane U) which shows a very regular dependence and a smooth anti-correlation 
with U until the collisional quenching of Clll]A1909 sets on. However, the IP of Si 
is just 8 eV. The majority of the A1814 doublet is emitted in the partially 



ionized zone (PIZ) near the hydrogen ionization front. Fig 10 (bottom 
panel) shows the behavior of the volume emissivity e times the geometrical depth 
within an emitting slab of g clS cLS du function of the depth itself. Since the actual line 
emission is proportional to e • h reported in Fig. [10} the emission of the AI111AI86O, 
and Sim] A1892 lines is negligible in the PIZ, but the one of H/3 certainly is not, 
and may be well dominating if A c ^> 10 23 cm~ 2 . Significant emission in the PIZ is 
expected also for the SillA1814 doublet. However, for this latter line, total emission 
should not depend strongly on column density, since the emissivity becomes very 
low at A c > 10 23 cm~ 2 . If higher N c are considered (up to 10 25 cm -2 ), ratios including 
Sim] A1892, AI111AI86O, and SillA1814 show a very weak dependence on column 
density, with changes a few percent at worst. 

We attempted to isolate a CivA1549 and a SilvA1397 component that 
corresponds to the AI111AI86O and Sim] A1892 lines. This can be a small part of 
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the total CivA1549 and of SilvA1397 emission, but there is no point to consider 
the whole CivA1549 emission when CivA1549 shows a large blueshift and is much 
broader than H/3 (Fig. 5 of Sulentic et al. 2000) and Sim] A1892 and AI111AI86O 
(Fig. 2 of Marziani et al. 2010). Taking into account various sources of ambiguity 
(mainly uncertainty in the quasar rest frame, S/N, blending with HeIIA1640), the 
CivA1549 BC component we measure with our fits is constrained within ± 50% 
at worst. Thus we find significant CivA1549 emission from the low-ionization gas, 
and the basic assumption is that its CivA1549 profile is the same as the AI111AI86O, 
Sim] A1892 lines. If one considers the emissivity behavior, CivA1549 and SilvA1397 
are obviously favored within the fully ionized zone. 

The ratio of CivA1549 and SilvA1397 over Cm]A1909 or over 
Silll]A1892 increases with ionization parameter in a way that is roughly 
independent on density until the collisional quenching of the semi- 
forbidden lines sets on (Fig. [9j). The previous considerations are helpful 
to understand why the CivA1549 and SilvA1397 ratios provide clear 
diagnostics of the ionization parameter and seem to be in (at least 
qualitative) agreement with ratios employing lower ionization lines like 
SillA1814/ Silll]A1892 that is also mainly sensitive to U. The striking 
fact that the CivA1549 and SilvA1397 emission confirms low-ionization 
supports the hypothesis that the four lines are emitted by the same 
region. The ratios involving Si only have the obvious advantage that 
the the determination of the physical parameters is not dependent on 
metallicity; the drawback of SilvA1397 is that this line is weaker and 
blended with Oiv]A1402. However, as already pointed out, the SilvA1397 
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BC intensity should be slightly affected by Oiv]A1402. 

7. Results on two Extreme, Elucidating Cases 

We analyze two extreme objects - one extreme Pop. A and one extreme Pop. 
B - with the same methodology we used previously. The aim of this analysis is to 
help the interpretation of the line components measured on the spectra of the 8 
z ~3 quasars. 

7.1. SDSS Weak Cin]A1909 source: SDSS J120144.36+011611.6 

Our project stems from the careful analysis of the I Zw 1 UV spectrum by 
Laor et al. (1997b), and the evidence of a well-defined trend in the A1900A blend 
along the El sequence (Aoki & Yoshida 1999, Wills et al. 1999, Bachev et al. 2004, 
Marziani et al., 2010). The NLSyl I Zw 1 is known to be a sort of extremum in the 
El sequence: it shows strong Fell and Fein, prominent AI111AI86O and Sim] A1892 
emission. It is an example of the A3 spectral type, whose median 1900A blend in 
shown in Fig. 3 of Bachev et al. (2004). Definitely, these sources are present also 
at intermediate to high redshift. We describe here the analysis of one source, SDSS 
J120144. 36+011611. 6 (Fig. [8]), which seems to be a high-redshift, high-luminosity 
replica of I Zw 1, with broader lines (the FWHM limit of NLSyls and Pop. A 
sources is luminosity dependent; see Netzer & Trahktenbrot 2007 and Marziani et 
al. 2009). These sources are Pop. A, and we assume that the profile of the BC is 
Lorentzian. Hereafter unshifted Lorentzian part of the line is said to be BC. Pop. 
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A sources are also free of any VBC, making the analysis of the blend simpler. 
Several considerations can be made from Fig. |HJ 

AI111AI86O The lines less ambiguous to measure are the AI111AI86O doublet 
because the two lines are less blended with other features, and they are remarkably 
strong. Sim] A1892 is more heavily blended with Clll]A1909 and Fein. However, 
the specfit routine allows usually a plausible deconvolution of Sim] A1892, making 
the AI111AI86O/ Sim] A1892 ratio very reliable. Fell and Fein are obviously strong. 
We can use Fell UV 191 to set roughly the level of Fell, while the feature at 2080A 
is helpful to estimate intensity of Fein. 

Fein Our specfit analysis produces a very weak Clll]A1909 component. A 
precise measurement of Clll]A1909 is cumbersome, since its intensity depends on 
the actual Fell I emission. The strong feature at A1914 could be dominating, and 
the template may seriously underestimate it (see Vestergaard and Wilkes 2001 
for several alternatives in the empirical Fein emission of I Zw 1). In any case 
the residual Clll]A1909 emission is small, suggesting that the A1900 blend, once 
believed to be mostly Clll]A1909, is actually almost void of Clll]A1909 emission in 
these A3-type sources. 

SillA1814 The SillA1814 line is well visible in the spectrum of SDSS 
J120144. 36+011611. 6 and can be used as a substitute of CivA1549 to mea- 
sure the ionization parameter. The ratio SillA1814/ Sim] A1892 is mainly sensitive 
to the ionization parameters, as it is the ratio CivA1549/ Sim] A1892. The ratio 
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SiiiA1814/ Sim] A1892 has the considerable advantage of being weekly dependent on 



metallicity. If metallicity is known (see also £8.1), the ratio ClvA1549/AllllA1860 
should be in principle preferred since AI111AI86O is emitted through a permitted 
resonance transition while SiiiA1814 is emitted in the collisionally excited, 
partially-ionized zone (PIZ). 



ClvA1549 We fit CivA1549 with a Lorentzian component that is unshifted + a 
CIVA1549 BBC. Note that the lack of shift in AI111AI86O and Sim] A1892 imposes 
a strong, determinant condition on the strength of the Lorentzian-component in 
CivA1549. It is important to stress that some CivA1549 BC emission is expected 
to be present according to our array of simulations. The ratio CivA1549/A1iiiA1860 
is rather modest, as we can easily see even by eye. We note in passing that the 
with of BCC is ~ 10000 km s" 1 , with a peak blueshift indicatively of —6000km s _1 , 
significantly larger than the one measured on the CivA1549 of I Zw 1 by Marziani 
et al. (2009). 



SilvA1397 We fit SilvA1397 with a Lorentzian component that is 
unshifted + a BBC (that may include SilvA1397 and Oiv]A1402 
contribution). Results that are consistent with the ones of CivA1549. 
The ratio SilvA1397/ Sim]A1892 has the advantage that is almost 
independent on metallicity. In principle, the crossing point between the 
SillA1814/ Silll]A1892 and SillA1814/ Silll]A1892 can set a metallicity 
independent point in the (n H , U) plane. If the accuracy of the SillA1814 
intensity is good then this point can be used to retrieve information on 
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metallicity (pTIj). 



We consider the measured line ratios in the plane U vs. rin (initally assuming 
metallicity equal to solar), and see where they cross. We find very high density and 



low ionization (Fig. 11) 



7.1.1. Along the El Sequence: A More Complex Scenario 



Looking at the fit solution (Fig. p[b)), it seems that our spectrum has almost 
no Cin]A1909. In many ways this is not surprising. The physical solutions in the (U, 



hh) plane of Fig. 11, points toward very low ionization (U ~ 10~ 3 — 1CT 2 ) and high 
density (nn> 10 12 cuT 3 ). At such high values of wh we expect that Cin]A1909 is 
collisionally quenched, and no significant emission. The ratio Ciii]A1909/CivA1549 
is expected to be just ~ 10~ 2 in the dense LIL-BLR. Therefore we can say that 
sources like SDSS and I Zw 1 are extreme because all of the A1900A blend is 
emitted by very low ionization, dense gas. 

However, as soon as we move away from spectral types A3+ along the El 
sequence, we see that the prominence of AI111AI86O diminishes greatly. The 
emission disappears altogether at the other end of the El sequence, where several 
lobe-dominated radio-loud sources are found. For most quasars we see that 
Cin]A1909 is rather strong and unmistakably present. At a first glance this 
complicates the interpretation of the spectrum. However, one has to consider that 
Clll]A1909 can be emitted only by gas of much lower % than that of the region 
where the bulk of Siiul A1892 and AI111AI86O is emitted. Our simulations show 
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that AllllA1860 intensity grows smoothly as a function of density in the ionization 
parameter range of interest — 3 < \ogU < — 1. This said, and considered the smooth 
trend we see from A3 to Bl + , the most reasonable conclusion is that a dense region 
emits significantly whenever AI111AI86O emission is detected, although the relative 
prominence of the dense BC changes along the El sequence: it accounts for the 
entire LIL emission only in the most extreme Pop. A sources. The sequence of 
Fig. 3 of Bachev et al. (2004) seems to be mainly a sequence of prominence of 
AI111AI86O+ Sim]A1892 vs Cm]A1909. 

So, a very important conclusion is that a very dense, low-ionization region 
exists in the wide majority of quasars. It is associated with Fell prominence, as 
such gas is expected to emit a strong low-ionization spectrum. The most extreme 
Fell emitters are also the most extreme AI111AI86O emitters; in some cases where no 
AI111AI86O emission is measured, we also fail to detect any Fell emission (Marziani 
et al. 2010). Where is this region located? Why there are such distinctive line 
profiles as Lorentzian? The second issue goes beyond the present paper; for the 
moment our aim is to measure line components of AI111AI86O, SillA1814, Sim] A1892 
and CivA1549 that come from this region and to estimate its distance from the 
central black hole (Q. 




3C 390.3 is a lobe-dominated (LD) RL source, with very broad emission lines 
and no optical Fell within detection limits (see Marziani et al. 2003, for the criterion 
used). It also shows a large peak shift in its H/3 profile, and prominent narrow lines 




7.2. 



The Other El Extremum: 3C 390.3 
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([Olll]A5007, H/3), but also narrow components of CivA1549 and Clll]A1909 which 
are all well separated from the broad H/3 profile. 

To deconvolve the blends around Clll]A1909 and CivA1549, we assume the 
same emission line profiles as for H/3. The H/3 broad profile can be described as the 
sum of a BC and a VBC. An unusual property of the BC is its large peak blueshift 
(it is unusual because of the shift amplitude: even if median spectra of Pop. B 
sources show a small BC redshift, there is a pretty large scatter with both red- and 
blueshifted peaks observed in individual sources; Zamfir et al. 2010). Applying in a 
self consistent way the BC and VBC to the Clll]A1909 and CivA1549 blends, leads 
to very interesting results as shown in Figure |8j 

• AI111AI86O is very weak or even absent within the S/N limits; 

• the ClvA1549 profile is very similar to the one of H/3: BC+VBC accounts 
for more than 90% emission, with a possible, very weak contribution of the 
blueshifted component which is usually dominating in Pop. A sources; 

• the red wing observed in the 1900A blend is accounted for only if there is a 
strong Clll]A1909 VBC, Fell and Fein emission being below the detectability 
limit in this source; 

• there is no evidence of a VBC in AI111AI86O; 

• the ratio CivA1549/Ciii]A1909 BC is ~ 10, a far cry from the ratios observed 
in extreme Pop. A sources like I Zw 1 and SDSS J120144.36+011611.6; 
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• even more interesting we find a CivA1549/Ciii]A1909 VBC «10, as for the 
BC. 

We conclude that this object is fundamentally different from Pop. A sources. 
The low ionization, high density BC seems to be absent. The similarity in the 
CivA1549/Ciii]A1909 ratios suggest that we are observing a gas in conditions similar 
to the one associated to the VBC. We predict that Fell will remain undetected or 
found to be weak even with S/N — > oo. 



In the plane (U, 71h) of Figure 11, the line ratios converge to a point at \ogU ~ 
-1.5, and log 72h~ 10.1 (the SillA1814 line cannot be measured accurately since 
S/N is poor). In this case, the Ciii]A1909/CivA1549 ratio also converges toward 
a (U, nn) value consistent with the one indicated by the AllllA1860/SillA1814. 
In most other cases, this does not happen because the AI111AI86O doublet is too 
strong to be accounted for by gas of 10 10 cm -3 even at super-solar metallicity. 
Emission occurs at pretty high ionization, in conditions that once upon a time were 
thought to be standard in quasars (Davidson and Netzer 1979). This Pop. A and B 
difference at the extrema was already pointed out, in a semi-quantitative way, by 
Sulentic et al. (2000) and Marziani et al. (2001). 

It is however important not to generalize the case of 3C 390.3 to the remaining 
Pop. B sources. In most of them, AI111AI86O is detected, and there is evidence 
of strong Clll]A1909. Actually Clll]A1909 emission seems to be appreciable in 
all of our VLT quasars. This means that we are in a composite situation, where 
low-ionization, high-density emission is present, along with significant VBC and 
other emission. The two extreme cases help us however to understand these more 
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complex cases. 



7.3. The Contribution of Lower Density Gas 

Once the true intensity of the BC components is retrieved, the 
presence of significant Clll]A1909 emission complicates the analysis. As 
pointed out, the photoionization solution for the BC suggests very high 
density, and in this region no Clll]A1909 emission is expected. Whenever 
Clll]A1909 is observed, we need to reverse the question: how much does 
any Clll]A1909 emitting gas contributes to the lines used for diagnostic 
ratios? Negligible contribution is expected to AI111AI86O. However, this 
is not true for CivA1549 and Silll]A1892. Especially among Pop. A2 
and A3 objects, it is not so obvious that the profile of Clll]A1909 and 
Silll]A1892 is the same. It could be well that the Clll]A1909 profile is 
narrower than the ones of Silll]A1892 and AI111AI86O (as found for SDSS 
J1201+0116), justifying the idea of Clll]A1909 emission from a disjoint 
region ( |6.2 ). For the BAL QSOs in our sample and AllllA1860-strong 



sources most of what we ascribe to Clll]A1909 could be actually Fein, as 
suggested by Hartig and Baldwin (1986). 

However, objects like the typical Al sources show Silll]A1892 and 
Clll]A1909 with similar profiles, so in the following we will assume the 
worst condition, that is all Clll]A1909 emitting gas is contributing to the 
BC lines. 
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To estimate the contribution of the Clll]A1909 emitting gas to 
Silll]A1892 and CivA1549, we can consider the trends observed along 
El. Silll]A1892 is strong when AI111AI86O is strong, and the Silll]A1892/ 
Clll]A1909 ratio is lower when Clll]A1909 is strong, as visible in both 
Baldwin et al. (1996) and Bachev et al. (2004), as well as in the present 
paper. The observed Sim] A1892/Clll] A1909 ratio in the Al median bin 
seems to be as low as in Mark 335, ~ 0.4. Pop. B sources that often 
show prominent Clll]A1909 also show a low Sim] A1892/Clll] A1909 ratio, 
as also appreciable in the median spectra of Bachev et al. (2004). These 
trends suggest that most Silll]A1892 is emitted where AI111AI86O is also 
emitted. As a consequence, any correction due to gas in different physical 
conditions (lower density) emitting Silll]A1892 is expected not to be 
dominant unless Clll]A1909 is extremely strong. The Clll]A1909 emitting 
gas should be at density logn ~ 9, or lower. Higher density would 
imply increasing Sim] A1892/Clll] A1909 to values that would exceed 
the observed ones even if the LIL-BC is not emitting any significant 



Sim]A1892 (see Fig. 13) 



7.3.1. Preliminary analysis of a low-z sample 

To set these trends on a more quantitative basis we considered 
the set of pre-coSTAR recalibrated sources, for which Clll]A1909 and 
CivA1549 data are publicly available (Kuraszkiewicz et al 2002, Evans 
& Koratkar 2004) .We performed measurements (interactively with the 
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task splot of iraf) for about 30 sources with the highest S/N, holding a 
Clll]A1909 blend that could be relatively easily deblended, and following 
the expectations described in § |4.4[ The rest frame equivalent width of 



AI111AI86O and Sim]A1892 are found to be highly correlated (Fig. [12] ) . 
The correlation is due to Silll]A1892 being stronger when AI111AI86O 
is strong, not necessarily because Clll]A1909 is strong: actually, the 
Silll]A1892/Clll]A1909 ratio can be low when Clll]A1909 is strong, and the 
AI111AI86O/ Silll]A1892 ratio achieves maximum values when Clll]A1909 is 
faintest. This is consistent with the bulk of AI111AI86O and Silll]A1892 
originating in the same region. At the same time the presence of 
Cm]A1909 lowers the AI111AI86O/ Silll]A1892, CivA1549/ Silll]A1892 and 
SilvA1397/ Silll]A1892 because of "excess" Silll]A1892 emission associated 
to the gas emitting Clll]A1909. 



7.3.2. The effect of low-density gas on the product (Uu-r) 

We expect that any correction will increase density (increase 
AI111AI86O/ Silll]A1892 lowering Silll]A1892) and decrease ionization 
parameter (lowering CivA1549 more than Sim] A1892), but that their 
product will be less affected. To show the amplitude of the effect we 
performed an experiment, adding to a pure, high-density solution, 



contribution from moderate density gas (log wh~ 9 — 10). The Fig. 14 
shows the displacement in the density — ionization parameter plane 
for several values of the Cm]A1909/ Silll]A1892 intensity ratio. The 
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values refer to the Clll]A1909 addition over Silll]A1892 in the ideal case, 
corresponding to observed ratios of 0.34, 0.7, 1.0, 1.1 for an addition of 
a Clll]A1909 component whose intensity is 0.4, 1.0 ,1.5 ,2.0 the intensity 
of the Silll]A1892 component associate to the high density gas. As one 
can see from the figure, even if deviation for U and % taken separately 
are significant, deviations for the product are by far less important. 
The largest change for the product is found to be ps 0.15, if we exclude 
the gray dots corresponding to log %= 10 (an unlikley case, since this 
would imply a correction Silll]A1892/Clll]A1909 > 1, inconsistent with 
what we observe when Clll]A1909 is strong). Following the expected line 



ratios of Fig. 13 we apply a correction to the BC fluxes that is 0.4 and 



1.5 Clll]A1909 for Silll]A1892 and CivA1549 respectively (Figures 15, 16 



and 11), corresponding to log £7 = —2 and lognH= 9. We remark that that 
if U is lower, the correction will have negligible effect, while assuming a 
larger U will lead to CivA1549 flux in excess to the one observed. 

To constrain the ionization parameter we can first consider that, 
since the Clll]A1909 gas comes from a (relatively) low density region, 
the contribution to SillA1814 is small: for logn ~ 9 and \ogU ~ —2, 
the contribution should be m 0.03 Clll]A1909. Second, another powerful 
feature is the SilvA1397 doublet (Baldwin et al. 1996): the SilvA1397 
doublet is less affected by the Clll]A1909 correction, the contribution 
from lower density gas being estimated w 0.25 Clll]A1909. The line ratio 
SilvA1397/AllllA1860 is also sensitive to ionization and less affected by 
any lower density correction (provided that the relative abundance of S 
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and Al stays the same, as it seems to be the case). 

In summary, corrected BC line intensities are computed as 
follows: / c (SillA1814) BC = J(SillA1814) B c - 0.03 I (Cm] A1909) BC ; ^ c ( 
Siin]A1892) B c = /( Sini]A1892) BC - 0.4 /(Ciii]A1909) bc ; J c (SiivA1397) BC = 
7(SiivA1397) BC - 0.26 J(Ciii]A1909) bc ; J c (CivA1549) BC = J(CivA1549) BC - 
1.5 J(Ciii]A1909) bc . 



8. Results on the z ~ 3 Quasars 

To estimate logn B and \ogU values, we use the cloudy contour plots of 
the ratios AI111AI86O/ Sim]A1892, SillA1814/ Silll]Al892, CivA1549/ Silll]A1892, 
SilvA1397/ Sim] A1892 showed in Fig. £^ The data points of our objects are 



in regions were the ratios are well-defined. The ratios CivA1549/ Sim] A1892, 
SilvA1397/ Silll]A1892, and SillA1814/ Sim] A1892 are mainly sensitive to the 
ionization parameter U, while AI111AI86O/ Sim] A1892 and Clll]A1909/ Silll]A1892 
are mainly sensitive to the electron density. We know that Clll]A1909 is collisionally 
quenched at log n e > 10 and in the contour plot for Clll]A1909/ Sim] A1892 we see a 
step around this value. 

We measure the BC intensity of Sim]A1892, AI111AI86O, SillA1814, SilvA1397 
and ClvA1549; with them we compute the diagnostic ratios (for 3C390.3 we use 



2 Note that there are regions where the ratio values are actually undefined: close to the high U 
limit (\ogU > -0.3), ratios AI111AI86O/ Sim]A1892 and SinA1814/ Sim] A1892(with n H < 10 9 ciir 3 ) 
should not be considered. 
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Cm] A1909 in place of SiiiA1814; however there is no object similar to 3C390.3 
among the z ~ 3 quasars). We present the fluxes of the line components in Tab. [3] 
and the equivalent width in Tab. [5] Table [4] shows the weak lines around CivA1549. 
For CivA1549 line we show the core, blue shifted and the very broad components. 
Errors are at a 2a confidence level, and include the sources of uncertainty described 



in £4J5 Errors are then quadratically propagated according to standard practice to 
compute intensity ratios and their logarithm. 

From Table [3] we can derive the diagnostic ratios. As we see in Fig. [TJ 
SillA1814 is absorbed by telluric B band in J00103-0037, J03036-0023, J20497-0554 
(most affected). We will not consider SiiiA1814 to compute un and U on those 
cases. However, if we take at face value the SillA1814 measure on J00103-0037 and 
J20497-0554, it will converge close to the point set by the remaining two ratios. 
We display on a graph a line representing the behavior of each ratio under the 
assumption of solar metallicity; the ideal point where the lines representing different 



diagnostic ratios cross determines the values of lognn and log U. Figs. [15] and [16 
shows the contour plots were we can see that the diagnostic ratios converge to rather 
well defined values. The cross point is very precise for the objects J00103-0037, 
J00521-1108, J02287+0002 (using z C iii), and J20497-0554 ; for the remaining 
objects J01225+1339, J02287+0002 (using z i), J02390-0038, J03036-0023, and 
J23509-0052 the cross point is slightly different. We must not forget the errors 
involving the fits, such as the changing of the shape profile that makes large the 
peak intensity if is Lorentzian or it could be less intense if it is Gaussian; the Fell 
pseudo-continuum contribution that affects principally to SiiiA1814 or the Fein that 
in some cases affects Clll]A1909. 
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In principle, the crossing point of the ratios SilvA1397/ Silll]A1892 
and SiiiA1814/ Silll]A1892 is independent on metallicity. Therefore, 
any significant disagreement between this crossing point and the ratios 
based on CivA1549 may indicate chemical composition different from the 
assumed solar one (§ |8.1 ). The difficulty here is the large uncertainty of 



the SillA1814 line. In all contour plots we show the ±2cr interval as a 
shaded band. So, it is proper to consider deviations from metallicity only 
in the case of 4 sources where the crossing point excluding SillA1814 is 
outside the uncertainty band. 

For the objects J00103-0037, J03036-0023 and J20497-0554 we exclude the 
SiiiA1814 line from the diagnostic ratios because it is affected by absorption. 
For the remaining objects, we take the average of the crossing contour plots 
of log(AllllA1860/ Sim] A1892) crossing with log(ClvA1549/ Sim] A1892) and 
log(AllllA1860/ Silll]A1892) crossing with log(SiiiA1814/ Sim] A1892). For 3C390.3 
we use log(Clll]A1909/ Silll]A1892) instead of log(SiiiA1814/ Sim] A1892). Table 
[6] summarizes the log % and log U values including their uncertainty. Since U 
and nn are not independent quantities (their correlation coefficient is found to 
be 0.55), we adopt the appropriate formula for the errors on the product Un-R 
(following Bevington 1969). We present average values of the crossing points for 



extreme objects of Fig. 11 and 1 18(d). In the SDSS J12014+0116 case we give full 
weight to the SillA1814 measurement. The crossing points disagree somewhat for 
SDSS J12014+0116. This could be due to an underestimate of both SillA1814 and 
CivA1549 in the fits. On the one hand, SillA1814 is clearly seen and strong but is 
contaminated by Fell blend; on the other hand, CivA1549 is strongly affected by 
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the blue-shifted component and by the assumption that it is of Gaussian shape. 
An increase by 30% in the measurement of the intensity of the two line would 
lead to a better agreement but this is somewhat an ad-hoc speculation. Rather, 
the significant difference between the crossing point of the SillA1814/ Sim] A1892 
and SilvA1397/ Sim] A1892 ratios and the other ones point toward strong metal 



enrichment. We will show in £8.1 that this is probably the case. 



At any rate, the convergence is toward a value of log U ~ — 3, lower than for 
the other z ~3 quasars. This is reflected in the CivA1549 EW of this source, also 
significantly lower. It is intriguing to note that the correction because of 
lower density drives the other sources toward values of U and rin that are 
closer to the ones of SDSS J12014+0116. 

Table [7] reports the values of the tblr and the Mbh of our 8 objects and the 
extreme objects in the last two rows. Column 1 identifies the quasar name; Col. 
2 gives the quasar proper distance in mega parsecs [Mpc]; Cols. 3 and 4 are the 
continuum specific flux value at 1350A and 1700A respectively, Col. 5 reports the 
FWHM in km s _1 for the broad components, Col. 6 is the Population designation. 
Cols. (7) to (10) report the logarithm of the size of the BLR in cm 
obtained from: a) 1Z Q , b) 1Z & line ratios corrected because of low density 
emission, c) 5Z , and d) 5Z Q line ratios corrected because of low density 
emission. Cols. (11) -(14) list the logarithm of the black hole mass in solar masses 
in the same order as for tblr- Finally Col. (15) is Mbh computed from Vestergaard 



and Peterson (2006) formula (Equation 10). We will explain in £10 how these 
quantities are computed. 
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8.1. Effects of Metallicity 

The strength of NvA1240 relative to CivA1549 and HellA1640 suggests 
supersolar chemical abundances (Hamann and Ferland 1993; Hamann & Ferland 
1999). Chemical abundances may be well 5 to 10 times solar (Dhanda et al. 2007), 
with Z ~ 5Z G reputed typical of high z quasars (Ferland et al. 1996). The El 
sequence seems to be mainly a sequence of ionization in the sense of a steady 
decrease in prominence of the low-ionization BC toward Population B (Marziani et 
al. 2001, 2010). However, this is not to neglect that metal-enrichment also plays a 
role, especially for the most extreme Pop. A sources i.e., those in bin A3 and higher 
(Sulentic et al. 2001). 

The lines employed in the present study come from carbon, silicon and 
aluminium; all these element can be significantly depleted from gas if dust grains 
are formed (e. g., Mathis 1990). However, the emitting regions where our lines are 
produced are thought too hot to contain significant amount of dust (a definition 
of broad line region is right the central engine region below the dust sublimation 
region: e.g., Elitzur 2009). In addition Si and Al are expected to be produced under 
similar circumstances in the late stage of evolution of massive stars (Clayton 1983, 
Ch. 7). We note also that the CivA1549/ Silll]A1892 and CivA1549/A1iiiA1860 
usually give results that are in perfect agreement in the plane {n^U). These 
findings support our assumption that, if metallicity variations are present, the 
relative abundance Al to Si remains constant. We considered two cases for 
enhanced metallicity: (l) constant solar abundance ratio Al:Si:C with 
Z = 5Z Q (5Z) ; (2) an overabundance of Si with respect to carbon by 
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a factor 3, again with Z = 5Z Q (5ZSA). This condition comes from the 
yields listed by Woosley and Weaver (1995) from type II Supernovae. 
The Si overabundance is also supported by the chemical composition of 
the gas returned to the interstellar medium by an evolved population 
with a top-loaded initial mass function simulated using Starburst 99 
(Leitherer et al. 1999). The abundance of Al should scale roughly with 
the one of Si. While some cases with Al scaling with C are possible from 
the Woosely and Weaver (1995) yields, they are rarer than cases in which 
Al scales with Si. This latter case is appropriate for the most massive 
progenitors. Also, the assumption of Al scaling with C with [Si/C] = 
0.477 would yield to implausible high density and lack of convergence 
to a well-defined solution for lognn ~ 14. We therefore assume in the 
following that Al scales with Si in the two cases listed above. 

An array of simulations as a function of ionization parameter and density was 
computed assuming the conditions (1) and (2) listed in the previous paragraph. As 
expected, if the solar metallicity is simply scaled by a factor (5Z) we find that the 
ratio AllIlA1860/ Sim] A1892 is not strongly dependent on Z: the ratio increases by 
about 40% passing from Z = 1Z & to Z = 5Z Q , for lognn~ 12 and logU «-2. The 
same is true for the SillA1814/ Sim]A1892 and SilvA1397/ Sim] A1892 ratios. Since 
the first ratio sets nu, and the last two U, the ratios mentioned in this paragraph 
should be preferred because they provide nn and U values that are weakly affected 
by a factor 5 change in metallicity. A posteriori we confirm that the effect on the 
product Un u derived also with ratios involving CivA1549 is negligible in case 5Z 
(and it should be even more so if a metallicity increase is Z Q < Z < 5Z Q ). 
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The two extreme cases seem to be revealing also as far as metallicity is 
concerned. Ratios converge to a fairly well defined point in the case of 3C 390.3 



(see Fig. 11, left panel, Z = 1Z & assumed). In this case there is no major 
evidence of supersolar metallicity. The converse is true in the case of SDSS 
J120144. 36+011611. 6. The ratios involving CivA1549 indicate a lower ionization 
level, with \ogU ~ —3. This is because C 3+ changes ionization state to C 2+ 
for smaller ionization parameters, and so CivA1549 rapidly disappears. 
The CivA1549 intensity depends weakly on Z while the AI111AI86O and Sim] A1892 
lines are more sensitive. The 5Z case yields values closer to the ones obtained from 
the Si and Al line ratios, but not yet in concordance. If we pass to case 5ZSA 
with a factor 3 Si overabundance, then the concordance of the line ratios is good, 
especially if no correction because of low density emission is applied. As a further 
confirmation we checked that the metallicity-dependent ClvA1549/SilvA1397 ratio 
is in very good agreement with the crossing point of the other lines. Therefore, 
in this case, we have independent measures of metallicity, tt-h, and U. The case 
of even higher metallicity, say Z ~ 1OZ , remains to be explored but may not be 
appropriate considering the good agreement with 5ZSA. 

In principle, the discrepancy in the intersection point, with Sim] A1892 / CivA1549 
and A1iiiA1860/CivA1549 yielding lower U than the Al-Si ratios, should signal 
a significant enrichment in Si and Al of the BLR gas. In this case the 
SilvA1397/ClvA1549 ratio should be helpful, as it can be assumed to be dependent 
mainly on the Si abundance relative to C. Here, more than precisely determining the 
exact abundance value we are interested in analyzing the effect of large metallicity 
changes on U, %, and their product. Appreciable discrepancy is visible in the 
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contour plots of the BAL quasar JO 1225+ 1339, J03036-0023, J20497-0554 and 
J23509-0052 if Z = Z Q is assumed. 

In the same plots made for Z = 5Z Q SA, the agreement becomes better (see Fig. 



18). High metallicity yields higher U and smaller % if emission line ratios involving 
ClvA1549 are considered. This reflect the increase in abudance of Si and Al with 
respect to C, and the fact the SillA1814, Sim] A1892, AI111AI86O lines are emitted 
at lower ionization than CivA1549. In the case of sources like J02390-0038 and 
J20497-0554, the discrepancy of the crossing lines might indicate Z & < Z < 5Z Q , 
more than the extreme enrichment like the one assumed in 5ZSA, as also suggested 
by the ClvA1549/SilvA1397 ratio. 



In Fig. 18 and 19 the contour plots are shown for the case 5ZSA. We consider 
uncorrected and corrected line ratios as two independent cases. The halftone bands 
also show the importance of accurate SiiiA1814 measurements to infer unambiguous 
constrain on metallicity (and, by extension, to estimate wh and U independently). 
The product n^U is however much less affected than wh and U individually (Tab. 

A first tblr and Mbh estimate can be obtained considering only the SiiiA1814/ 
Sim] A1892 and AI111AI86O/ Sim] A1892 ratios. We conclude that the effect of 
scaling the metallicity up to Z = 5Z Q is within the uncertainty of the method, 
and particularly small if the ratios SillA1814/ Silll]A1892, AI111AI86O/ Silll]A1892, 
and SilvA1397/ Sim] A1892 are considered to compute nn and U. It is significant 
if strong enrichment of Al and Si over C occurs. A more refined approach could 
exploit the dependence of AllllA1860/ClvA1549 (and Sim] A1892/CivA1549 and 
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SilvA1397/ClvA1549) on Z to build a 3D diagram where n H , U, Z along with Si-Al 
enrichment can be determined independently. 

8.2. Determining the best estimate of n H , U 

There are two sets of line ratios for each object: one comes from 
the specfit results, and the other is the one computed after correcting 
the specfit results because of low density contributions. Solutions 
with corrected value usually point toward very high density and 
low ionization, predicting SillA1814 emission even twice as strong as 
Silll]A1892. However, increasing the metallicity for corrected ratios leads 
to better agreement among lines and more reasonable tih, U values, 
while leaving the product fairly unaffected. Especially the assumption 
of Silicon - Aluminium enrichment improves the concordance in the 
intriguing cases of SDSS J12014+0116, J01225+1339, J02287+0002 
(using z oi ), the extreme objects and the two BAL QSOs in our sample. 
However, apart from the case of SDSS J12014+0116, the enrichment 
is probably excessive. From the discussion above the independent 
determination of % and U seems possible only if metallicity is at 
least roughly known. We exclude metallicity cases where we find a 
sizeable disagreement in the crossing points (with the exceptions of 
ratios involving SillA1814: the ionization level can be estimated in a 
Z-independent way using the SilvA1397/SillA1814 and the SilvA1397). 
We consider each individual source with line intensity before and after 
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correction to obtain two independent sets of product n^U values. As 
mentioned, changing metallicity is not affecting the product null as 
much as ?t,h and U individually. In case of concordance of the crossing 
points and of high accuracy in the SillA1814 ratio, n H U, and Z can be 
considered independently determined. In Table [6] we indicate the values 
that are deemed most appropriate. 



9. A Photoionization Method to Compute the Broad Line Region 
Distance and the Black Hole Mass. 

The distance of the broad line region (Vblr) from the central continuum source 
and the black hole mass (Mbh) are key parameters that let us understand the 
dynamics of the gas in the emitting region and the quasar behavior and evolution. 
In this work we will use a method based on the determination of and U to 
compute tblr- Eq. [T] can be rewritten as 



f + °° ^dv 



1/2 



and also as 



r BLR 



1/2 



d r , 



(3) 



(4) 



where h is the Plank constant, c is the light speed, d p is the proper distance. 
The integral is carried out from the Lyman limit to the shortest wavelengths on 
the rest frame specific flux f\. For the integral we will use two Spectral Energy 
Distributions (SEDs): one described by Mathews & Ferland (1987) and one by Laor 



et al. (1997a), also reproduced in Fig. 20 
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Eq. [4] becomes: 



r BLR « 93 ■ (C/n e ) 1Q 5 ■ fl 15 Qlo.iC(z, 0.3, 0.7) Id (5) 



where Qh = J^ hy n\Xd\, and £ (2, 0.3, 0.7) is an interpolation function 
for dp as a function of redshift. 

Qh is 0.00963 cmA in the case the continuum of Laor et al. (1997) 
is considered; Q H ^0.02181 cmA for Mathews & Ferland (1987). We use 
their average value, since the derived U and n-^ are not sensitive to the 
two different shapes to a first approximation]^] 

Knowing tblr we can calculate the Mbh assuming virial motions of the gas 

Mbh = f q • (6) 

or, 

M BH = ^fo.7 5 (FWHM) 2 r BLR (7) 

with the geometry term / w 0.75, corresponding to /0.75 ~ 1.0 (Graham 
et al. 2011, see also Onken et al. 2004 and Woo et al. 2010). Collin et 

al. (2006) suggest that / is significantly different for Pop. A and B sources; we do 
not consider here their important result for the sake of comparison with previous 
work ( ^10.2 ). The resulting tblr and Mbh are reported in Columns 7 to 14 of 



Table [7} Errors are at 2cr confidence level and have been computed propagating 



3 Since the Laor et al. (1997) continuum produces a fewer ionizing photons, the same value of 
U is obtained at a smaller distance. 
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quadratically the major sources of uncertainty. More precisely, in addition to 
the error on login^U), the tblr determination is affected by the uncertainty in 
the spectrophotometry (specific fluxes of Col. 3 of Tab. [7]), and errors on the 
shape of the ionizing continuum. The two SEDs that we assumed as extreme 
yield a difference in ionizing photons of a factor 2.2. At a 2er confidence level this 
corresponds to an uncertainty of ±0.065 in logarithm. An additional source of 
uncertainty affects Mbh due to the FWHM determination. Errors on FWHM are 
quadratically added to the uncertainty on t B lr in the values reported as logM B H 
errors. 

10. Discussion 

10.1. Previous work 

There have been several studies aimed at computing r blr and M BH . A direct 
measure of tblr through reverberation mapping requires an enormous amount 
of observational effort and has only been applied to a relatively small number 
of quasars: slightly less than 50 objects with z < 0.4 (Kaspi et al. 2000, 2005, 
Peterson et al. 2004; Bentz et al. 2010). A second way to measure tblr uses a less 
direct method. Kaspi et al. (2000, 2005) and Bentz et al. (2009) used reverberation 
mapping results to find, in an empirical way, a relationship between tblr and the 
optical continuum luminosity at 5100A, 

r BL R oc L a (8) 
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with a ~ 0.52. Vestergaard & Peterson (2006) obtained a similar result for the 
optical continuum luminosity with an a ~ 0.50 and for the UV continuum at 
1350A, a ~ 0.53. These relationships have been used to compute the tblr not only 
for nearby objects, but also for high redshift, high luminosity objects. There are 
other works that use single epoch spectra and the continuum at 3000A, obtaining 
an a « 0.47 (McLure & Jarvis 2002). 



We can rewrite Eq. [6] as 



FWHM 2 L ft 

M BH oc / . (9) 



H/3 is a low ionization strong line whose FWHM has been widely used 
to determine the Mbh for objects mainly up to z < 0.9; above this limit IR 
spectrometers and large telescopes are needed to cover the redshifted line. For 
distant objects (z ~ 2), an alternative is to use CivA1549, a high ionization line 
emitted in the UV. However, this line should be used with caution because the 
line is often blueshifted. This means that at least part of this line is likely emitted 
in an outflow (Sulentic et al. 2007; Richards et al. 2010). Thus the estimation of 
Mbh using FWHM (CivA1549) tend to be systematically higher than those using 
FWHM(H/3), especially for objects of Population A. 

10.2. Comparison with Vestergaard and Peterson (2006) 

Vestergaard & Peterson (2006) used the relationship tblr oc L - 53 to obtain the 
following formula that relates Mbh to the FWHM(CivA1549) and the continuum 
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luminosity at 1350 A: 

2 r \ r horn X \ "1 0-53" 



logM BH (Civ) = log 



FWHM(Civ) 
1000 km s" 1 



AL A (1350A) 
10 44 ergs s- 



■ + (6.66±0.01)-s f . (10) 



The scale factor s f m —0.27 sets the masses to the / value obtained 
by Graham et al. (2011). In Cols. 11 and 15 of Table [7] and in Fig. 



21 



we 



compare our Mbh results with those using Eq. [TOj We do not apply corrections 
for radiation-pressure effects that are likely relevant especially for objects radiating 
at large Eddington ratio (Netzer 2009; Netzer & Marziani 2010). The difference 
between this computation and the one reported in Sulentic et al. (2007) is that 
in the latter work the blueshifted component was not separated from the broad 
component of CivA1549 just to show how larger values of FWHM(CivA1549) 
yielded Mbh much larger than the ones derived from FWHM(H/3) in Pop. A 
objects. 

We compare the masses obtained using our photoionization method with those 



of Vestergaard & Peterson (2006) in Fig. 21 We use the FWHM of the BC as an 



estimator of the virial line broadening. From our results we can see that the 
masses agree within less than la uncertainty in the luminosity correlation 
(0.33). There is a systematic offset of 0.17±0.10 if uncorrected ratios 
are used. The mass values obtained after correction for low density gas 
are systematically lower. This happens because the correction increases 
the product UriR, lowering tblr and hence Mbh- The Mbh obtained after 
corrections are within the error bars. The systematic offset is then 
0.13 ± 0.12. It is important to consider that the computed correction is 
in many ways a maximum correction. Clll]A1909 emission is assumed 
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to have the same FWHM of Silll]A1892 while it could be significantly 
narrower; in addition part of the Clll]A1909 emission could be due to 
Fein A 1914. In many sources of Pop. A the correction could be ignored 
altogether. The present results indicate that the Kaspi et al. (2000) relationships 
can be extended to be used in high redshift objects (or at least until z~3) if the 
FWHM of the core broad line region can be well determined and measured. In 
order to do this, we need: 

• Spectra with S/N high enough to see the profile shape that allows 
decomposition of the CivA1549 line, especially to separate the blue component 
from the broad core; 

• to follow the methodological considerations explained in §[4] . 



Fig. 21 should be looked at with two cautions. First, the correlation is 
dominated by the luminosity dependence of tblr, used to compute Mbh in both 
cases. Second, the spread of M B h values is small, less than one order of magnitude 
(and most objects have statistically indistinguishable masses). Our estimated error 
bars are however smaller compared to the spread expected on the basis of the 
tblr-L correlation which is, according to Vestergaard & Peterson (2006), ±0.66 at a 



2a confidence level. The two shaded bands of Fig. 21 limit the region where we can 
expect to find data points on the basis of the tblr--^ correlation. Clearly, a proper 
interpretation of the CivA1549 profile may help to reduce the scatter. In any case, 
our method should provide Mbh estimates with somewhat lower uncertainty. It is 
interesting to note that SDSS J12014+0116 appear at the largest Mbh, and the 
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agreement with the r blr - -^ is very good. The r blr value of 3C 390.3 we obtain is 
fairly uncertain due to the low S/N. The size computed in the present paper is 
larger than the reverberation mapping derived r B LR for H/3 and CivA1549, 
although consistent with the value derived from Lya. 3C 390.3 has the 
unusual property of having a response time longer in CivA1549 than 
in H/3, although the large error bars do not exclude that the two lines 
respond with similar time. This behavior might be related to different 
physical conditions found for this object. 

10.3. The LIL-BLR 

While the existence of high density and low ionization has been invoked since 
long to explain Fell emission (especially by S. Collin and collaborators, as mentioned 
in the introduction), we have provided additional evidence that high density and low 
ionization are indeed diagnosed from emission lines other than Fell blends, and that 
conditions are the one producing most line emission in extreme Pop. A quasars. Last 
but not least, both U and wh can be observationally determined with reasonable 
accuracy from the diagnostic ratios SiiiA1814/ Sim]A1892, SilvA1397/ Silll]A1892, 
AhllA1860/ Sim]A1892, and CivA1549/(A1iiiA1860 or Sim] A1892). Metallicity can 
be also constrained from the previous ratios as well as from ClvA1549/SilvA1397. 
These ratios are pretty well defined, while estimating Clll]A1909 and FelllA1914 
relative contribution is not relevant to our method (with the exceptions of sources 
like 3C 390.3). This low ionization BLR (or LIL BLR) has very similar properties 
to the Ol and Call emitting region identified by Matsuoka et al. (2008). The 
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LIL-BLR seems to be present in the vast majority of quasars, probably all the ones 
with significant Fell emission (Marziani et al. 2010). The low values of U could be 
a consequence of the high density rather than of a far away location of the emitting 
region. The assumption of a single well defined value of U and uh is probably an 
idealization even for the LIL-BLR taken alone; however the convergence of emission 
line ratios toward a well defined point, along with the ability to qualitatively explain 
Fell and most H/3 emission, indicate that the LIL-BLR might be a region with a 
small range of low U and high %. 

Do density and ionization parameter in the LIL BLR converge to a single well 
defined value with a small dispersion? A tentative answer comes from the values 
reported in Tab. [6j Excluding 3C390.3, the range U and «h span is not very large 
(even considering changes in metallicity): less than one order of magnitude around 
log Wh~ 12.5, and ionization parameter ~ —2.75, and the spread is not much 
larger than the uncertainties in the individual measurements. The product U 
seems to be fairly stable. We have < log(C/-riH)>~ 9.5, with a sample dispersion 
of 0.15 (excluding 3C 390.3). We applied the same method to 14 low- z quasars 
(Negrete et al. 2010), and we obtain < log([/-rtH)>~ 9.7, with a dispersion of 
0.3. Assuming log(£/-nH) ~ 9.6 could be a good approach to estimate Mbh if 
elaborate measurements on CivA1549 and Cin]A1909 are not possible, and only 
a rough estimate of the Sim] A1892 BC FWHM is available. On the other hand, 
there is still a major effect of measurement errors on the uncertainty derived for 
log([/-7iH); instrumental improvements may lead to a significant appreciation of 
object-by-object diversity. 
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The U-uyi values reported in Table [6] are not very far from the average value 
obtained by Padovani and Rafanelli (1988). This is not surprising since the spectra 
at low and high z seems to show the same diversity, classified through the El 
sequence and the Pop. A/Pop. B distinction. In other words, NLSyl-like sources 
whose spectrum is similar to I Zwl appear to be present at high redshift, meaning 
that the LIL-BLR remains strong and prominent over a wide range of redshifts. A 
second reason is that Padovani and Rafanelli (1988) considered H/3. Emission of H/3 
can be significant under a much wider range of U and %; however, the AI111AI86O, 
MgllA2800, Fell emitting region should be also a strong producer of H/3. This region 
with well-defined physical conditions is expected to be the emitter of the core of 
H/3, i.e., the part of the line responding more strongly to continuum changes. 

10.3.1. Verification on EW and Line Luminosity 

Considering that we have very low ionization parameter values, a legitimate 
question is whether we have a sufficient number of photons to explain the EW and 
luminosity of the emission lines. We made a preliminary check of consistency for the 
EW from CLOUDY simulations. A second, a-posteriori test was to consider the 
predicted line luminosity assuming the actual luminosity of the quasar, the density 
and the distance tblr derived from our method, and spherical geometry. 

A remarkable property common to all spectra is the low EW of the emission 
lines. The extreme source SDSS J120144.36+011611.6 has total jy(CrvA1549) w 
19 A, which becomes ~ 7 A if only the LIL component is considered. The whole 
1900 A blend has W « 20 A; the individual EWs of Sim] A1892 and AI111AI86O are 
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just a few A. Similar considerations apply to the other high-;? quasars, especially 
for Pop. A sources. In this case we have W / (CivA1549)~ 20 — 30 A, the equivalent 
width of the whole 1900 A 30 — 40A including the uncertain contribution of Fein 
and Cin]A1909. This has the important implication that the ionization parameter 
cannot be very large. From our simulations, we deduce that the LIL-BLR log U is 
always < —2.0. On the other hand, toward the low U, high density limit emission 
lines tend to disappear altogether. The predicted EWs become too low to account 
for the observed EW with a covering factor f c < 0.5 if log £7 < —3. The values 
obtained after correction log £7 ~ —3.25 are still possible within the condition of 
fc J$ 0.5 since the CivA1549 EW is greatly diminished. Cases of very low U are 
rather indicative of strong metal enhancement than of extremely low ionization 



level (see £8.1). If we take log 71h~12.5, logU = —2.75, the predicted equivalent 
width of is W / (CivA1549)~37 A for a covering factor 0.5, close to the largest values 
observed in our quasars. However, The CivA1549 EW of SDSS J12014+0116 can 
be accounted for by the values of logU ~ —3.00 and log 71h~ 12.75 - 12.50. 

If we model the BLR with a spherical geometry where the emitting gas covers 
a fraction f c of the continuum, CLOUDY computations confirm that the line 
luminosity can be accounted for. The luminosity at 1700 A of J00 103-0037 is 
logAL A ~ 46.6; the predicted luminosity of CivA1549 is logL(ClvA1549) ps 44.5, 
under the assumption of Z = Z®, f c = 0.1, and Mathews & Ferland continuum 
shape. The observed CivA1549 line luminosity of J00103-0037, « 9 ■ 10 44 ergs s~\ 
is obtained with f c w 0.3. 



10.3.2. Analogy with r] CarincE 



The physical conditions we envisage for the BC of quasars find a correspondence 
in the so-called Weigelt blobs of 77 Carinae, located in the equatorial plane of 
the system, perpendicular to the symmetry axis of the bipolar lobes forming the 
homunculus nebula (cf. Marziani et al. 2010). Unlike the gas of the bipolar 
lobes, predominantly shock heated, the Weigelt blobs are believed to be dense 
gas photoionized by the radiation associated to the central, massive star and to a 
possible companion (e.g., Johansson et al. 2000; Davidson 2005). The spectrum 
of the Weigelt blobs shows very weak Clll]A1909 along with a prominent line 
at A1914, ascribed to the z 7 P® — > a 7 Ss Fein transition. The line appears very 
strong because the upper level is populated by Lya fluorescence. This very same 
process is expected to be present also in quasars. Indeed, in I Zw 1, where lines 
are narrow, and in SDSS J120144. 36+011611. 6 the peak emission at around 1914 
A is actually visible. The amount of Lya pumping to the upper level (z 7 P§) of 
the UV 34 cannot be estimated through the standard edition of CLOUDY (the 
relevant levels of the UV 34 multiplets of Fe +2 ion are not included). Additional 
photoionization computations including a suitable Fe +2 model and line transfer 
should be considered. This is beyond the aim of the present study; we can conclude 
in a qualitative fashion that the spectrum of the i] Carinae blobs supports a view of 
the low-ionization part of the BLR that is not conventional: very high density gas, 
at very low ionization. 



-67- 



11. Conclusions 

In this paper we presented new observations of eight high redshift quasars. 
The spectra were meant to provide high S/N, moderate resolution data on which 
the CIVA1549, Sim]A1892, AI111AI86O, and SiiiA1814 emission line profiles could 
be accurately analyzed. Line profile fits allowed us to isolate a specific component 
whose intensity ratios were used to derive consistent values for electron density 
and ionization parameter. This line component (LIL BC) seems to be emitted 
predominantly by low ionization, high density gas in the majority of quasars 
studied thus far by us. 

These results permitted us to compute the product tir-U and hence the size 
of the Broad Line Region and the central black hole mass. The method described 
in this paper rests on the assumption of photoionization as the mechanism of gas 
heating; on the assumption of isotropic luminosity, and on line ratios predicted by 
CLOUDY simulations. The photoionization method explored in this paper offers an 
estimate of r blr for each quasar, with some advantages on the r blr valued derived 
from the luminosity-size correlation. The luminosity correlation suffers from large 
scatter and is simply extrapolated to very high luminosity without any support since 
there are, unfortunately, no conclusive results on reverberation of high luminosity 
quasars even if heroic efforts are underway (e.g., Trevese et al. 2007, Botti et al. 
2010). We found that the black hole masses derived from the computed tblr and 
from the virial assumption are in good agreement with the ones derived from the 



luminosity-size relationship. Actually, Fig. 21 suggests that we might have reduced 



the errors of the Mbh by a factor of two with respect to the expectation from the 
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Vestergaard & Peterson (2006) relationship. 

We repeat that our Mbh and tblr results are based on the product n^-U and 
not on values of and of U taken separately. It seems that this product converges 
to two typical ranges of values, one of them associated to low-ionization, high 
density gas (the LIL-BLR). For our % and U determinations we do not use ratio 



Clll]A1909/ Silll]A1892 except for 3C 390.3. As we discussed in this ratio 
should not be considered at high density because Clll]A1909 is collisionally quenced 
if riH > 10 10 cm -3 . Cin]A1909 is produced in conditions that are very different 
from the ones we found for the LIL-BLR. While the method can be applied to most 
quasars, the application seems to be especially straightforward to quasars whose 
spectrum is like SDSS J120144. 36+011611. 6 (if high metallicity is properly taken 
into account) or, at the other end, 3C390.3. In the first case we have dominance by 
the LIL-BLR, in the second case the LIL-BLR seems to be completely absent and 
physical conditions look radically different: high ionization and moderate density. 
An inspection of SDSS spectra covering both the 1900 A blend and CivA1549 
(up to z ~ 3.5) shows that SDSS J120144. 36+011611. 6 has many replicas at high 
redshift, accounting for at least a few percent of all quasars. These high-metallicity 
objects should be the first candidates to expand black hole mass computations to 
high redshift without relying on the tblr - L correlation. 

To apply the photoionization method in the most effective way, determining nn 
and U with the lowest uncertainty, spectral data should be of moderate resolution 
(A/AA ~ 1000) as well as of high S/N. If the SillA1814 line can be measured in an 
accurate way, it would be possible to derive independent estimates of U, n H , and 
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Z/Z & in most quasars. 

Especially the most extreme (in terms of AI111AI86O strength) objects 
in bin A3 and A2 hold the promise to make possible an independent 
estimate of n H , U, and metallicity. Clearly, objects in bin Al resembling 
their median spectrum are not well suited for an application of the 
method. Also, any source with Cin]A1909/ Silll]A1892> 1 is subject to a 
large correction. In light of the many uncertainty, an average value of the 
product Un (obtained from the objects of the other spectral types) could 
be considered. Pop. B objects should not avoided entirely, especially 
whenever Silll]A1892 > Cm]A1909 after VBC removal. 

The present exploratory analysis emphasized several sources 
of uncertainty. However, the parameter needed for r B LR and M B h 
computation, the product Un H , seems to be fairly stable and well-defined. 
Even with an error of a 0.3 in logarithm, the square root will be subject 
to a 0.15 uncertainty in logarithm, much lower than the uncertainty 
associated with the rBLRluminosity correlation. The large intrinsic spread 
of the correlation at low luminosity, its uncertain extrapolation at very 
high luminosity make preferable a one-by-one determination based on 
physical properties of an emitting region that remains similar to itself. 
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Table 1. Line Components in the A 1900 blend. 



Ion 


A 


X 


Ei - 




Transition 




n 


c 


Note 




A 


eV 


eV 




S" 1 


cm 


-3 




Si II 


1808.00 


8.15 


0.000 - 


6.857 




2.54- 10 6 






1 


Si II 


1816.92 


8.15 


0.036 - 


6.859 




2.65 • 10 6 






1 


Al III 


1854.716 


18.83 


0.000 - 


6.685 


2p S/2 ^1/2 


5.40 • 10 8 






1 


Al III 


1862.790 


18.83 


0.000 - 


6.656 


2 po ,2c , 
M/2 *l/2 


5.33 • 10 8 






1 


[Si III] 


1882.7 


16.34 


0.000 - 


6.585 


3 P 2 ° -> 1 S 


0.012 


6.4 • 


10 4 


1,2,3 


Si III] 


1892.03 


16.34 


0.000 - 


6.553 


3 P° ->■ ^0 


16700 


2.1 • 


10 11 


1,4,5 


[C III] 


1906.7 


24.38 


0.000 - 


6.502 


3 P 2 ° ->■ 1 S 


0.0052 


7.7- 


10 4 


1,2,6 


C III] 


1908.734 


24.38 


0.000 - 


6.495 




114 


1.4- 


10 1() 


1,2,4,5 


Fe III 


1914.066 


16.18 


3.727 - 


10.200 


2 7 P 3 ° -> a 7 S 3 


6.6 • 10 s 






7 



Note. — All wavelengths are in vacuum. (1) Ralchcnko, Yu., Kramida, A.E., Reader, J., 
and NIST ASD Team (2008). NIST Atomic Spectra Database (version 3.1.5). Available at: 
http://physics.nist.gov/asd3. 2: Feibelman & Aller (1987). 3: n c computed following Shaw & 
Dufour (1995). 4: Morton (1991). 5: Fcldman (1992). 6: Zheng (1988). 7: Wavelength and A ki 
from Ekbcrg (1993), energy levels from Edlen and Swings (1942). 
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Table 2. Basic Properties of Sources and Log of Observations. 



Object name 


m B 


z 


Line 


M B 


Flux 6cm (mjy) 


Date 


DIT 


N exp 


Airmass 


S/N 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


J00103-0037 


18.39 


3.1546 


1 


-25.68 


0.40 


2006-11-08 


1139 


3 


1.16, 1.13, 1.11 


60 


J00521-1108 


18.70 


3.2364 


2 


-25.39 


0.43 


2007-01-01 


1199 


3 


1.15, 1.21, 1.29 


41 


J01225+1339 


18.24 


3.0511 


1 


-25.80 


* 


2006-11-08 


1259 


2 


1.36, 1.32 


92 


J02287+0002 


18.20 


2.7282 
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-25.72 


0.35 


2006-12-16 


1259 
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1.10, 1.12 


67 


J02390-0038 
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-25.36 
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57 
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88 


J20497-0554 
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-25.79 
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54 


J23509-0052 


18.67 


3.0305 
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-25.36 


0.41 
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1.10, 1.11, 1.14 


62 



*Not in FIRST 
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Table 4. Weak lines around CivA1549. a 



Object 




Niv]A1486 


SiiiA1533 


HeoA1640 










BC 


blue 


J00103-0037 




2.4 ± 1.9 


1.1 ± 0.8 


1.3 ± 0.4 


2.0 ± 1.8 


J00521-1108 




0.1 ± 0.2 


1.1 ± 1.1 


1.6 ± 1.2 


0.1 ± 0.2 


J01225+1339 






1.0 ± 1.8 


3.2 ± 3.4 


4.3 ± 1.8 


J02287+0002 


(i) 




1.0 ± 0.8 


0.3 ± 0.1 


0.2 ± 0.2 


J02287+0002 


(2) 




0.9 ± 0.8 


0.2 ± 0.1 


0.0 ± 0.2 


J02390-0038 




0.5 ± 0.9 


0.9 ± 0.9 


0.2 ± 0.4 


1.6 ± 1.1 


J03036-0023 






1.5 ± 1.2 


0.5 ± 0.5 


9.2 ± 3.4 


J20497-0554 




0.3 ± 0.6 


1.5 ± 1.5 


1.8 ± 0.8 


4.6 ± 2.7 


J23509-0052 






0.4 ± 0.4 


0.2 ± 0.9 


2.5 ± 1.6 



Extreme Objects 

J12014+01161 3.0 ± 3.5 4.9 ± 5.8 1.5 ± 1.7 10.3 ± 9.2 

3C 390.3 5.6 ± 2.7 1.6 ± 1.5 2.8 ±1.5 



Note. — (a) Units arc 10~ 14 ergs s _1 cm -2 A -1 . (1) Considering 
z Ol\l'Mi- (2) Considering 2C-T//]A1909- We do not show HcllA1640yBC' 
because is very weak, when is considered. 
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Fig. 1. — Sample of 8 VLT spectra. Abscissa is obseved wavelength in A, ordinate 
is specific flux in units 10~ 16 ergs s" 1 cm -2 A -1 corrected for Milky Way Galactic 
extinction. The superimposed dotted line is before atmospheric bands subtraction. 
We show the positions of the lines of our interest Cin]A1909, Sim] A1892, AI111AI86O, 
SillA1814, FellA1787, CivA1549 and SilvA1397. J01225+1339 and J02287 are BAL 
quasars. 
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Fig. 2. — OlA1304.8 used to place the restframe showed by the mark. Abscissa is 
rest frame wavelength in A, ordinate is specific flux in units 1CT 16 ergs s _1 cm -2 A -1 
corrected for Milky Way Galactic extinction. In J00521-1108 and J02390-0038 the 
peak of OlA1304.8 is not observed clearly. In J00103-0037, J02287+0002 and J20497- 
0554, the redshift results, using both OlA1304.8 and Clll]A1909, are ambiguous. In 
J01225+1339, J03036-0023 and J23509-0052, the redshifts obtained using OlA1304.8 
or Clll]A1909 are consistent. 




Fig. 3. — Sample of 8 VLT spectra in rest frame wavelength. Abscissa is rest frame 
in A, ordinate is specific flux in the rest frame in units 1CT 13 ergs s _1 cm~ 2 A -1 . 
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Fig. 4.— Fits for Pop. A objects: J03036-0023 (a, b), J20497-0554 (c, d), J23509- 
0052 (e, f). Upper panels show the fits and the lower pannels under the fits show the 
residuals and also the fitted absorptions lines. Upper abscissa is rest frame wave- 
length in A, lower abscissa is in velocity units, ordinate is specific flux in arbitrary 
units. Vertical dashed line is the restframe for CivA1549 and Clll]A1909. Long 
dashed line is the fit, solid dark lines are the broad components: CivA1549 in left 
panels and Clll]A1909, Silll]A1892, AI111AI86O, SillA1814 in right panels. Dotted 
dark lines under AI111AI86O show the doublet. Short dashed line is Fell. Fein is 
shown in dash-triple-dot line in the right panels. Dash-dot line in the left panels is 
the blue-shifted component of CivA1549 while dotted line is the very broad com- 
ponent, present also in Clll]A1909 for Pop. B objects. In the left panels we show 
with faint lines the contribution of NivA1486, SillA1533 and HellA1640 core and 
blue-shifted components. For colors see online figures. 
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Fig. 5.— Fits for Pop. B objects: up J00103-0037, middle J00521-1108, low J02390- 
0038. Units and meaning of symbols are the same of Fig. |4} 
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Fig. 6.— Fits for BAL quasars: up J01225+1339, middle J02287+0002 using z OiX130i 
rest frame, low J02287+0002 using zcmAi909 rest frame. Note in J02287+0002 the line 
displacement with the consequently line intensity changes, specially in Clll]A1909, 
Sim] A1892 and CivA1549 broad and blue-shifted components. Units and symbols 
are the same as in Fig. |4} 
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Fig. 7.— Fits for SilvA1397: (a) J03036-0023, (b) J20497-0554, (c) 23509-0052, 
(d) J00103-0037, (e) J02390-0038, (f) J01225+1339, (g) J02287+0002 and (h) 
J12014+0116. We do not measure SilvA1397 for J00521-1108 and 3C 390.3 because 
they have low S/N. We asume that the profile of SilvA1397 can be fitted with the 
same three possible components (core, blue-shifted and red-shifted) as in the case of 
CivA1549. Units and symbols are the same as in Fig. |4| 
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Fig. 8.— Fits for the two extreme objects: up J120144.36+011611.6, low 3C390.3. 
Units and meaning of symbols are the same of Fig. |4| For 3C390.3 it was needed 
to fit a narrow unshifted component that we show in dash-dot-dot line. Note for 
3C390.3 SillA1814 is almost absent. 
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Fig. 9.— Isocontours for the ratios (a) log(Cm] A1909/ Sim] A1892), (b) 
log(AlmA1860/ Sim]A1892), (c) log(SinA1814/ Sim] A1892), (d) log(SiivA1397/ 
Sim]A1892), (e) log(ClvA1549/AllllA1860), and (f) log(ClvA1549/ Sim]A1892) de- 
rived from CLOUDY simulations. Abscissa is electron density in cm -3 , ordinate is 
the ionization parameter, both in logarithm scale. 
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10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 11 10 12 
h [cm] 

Fig. 10. — Top: Ionic Fractions as a function of the geometric depth h in the gas 
slab. Bottom: Line emissivity per unit volume in units of ergs s -1 cm -3 multiplied 
by depth. 
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Fig. 11.— Contour plots for the extreme objects (a) 3C390.9 and (b) J12014+0116. 
Abscissa is electron density in cm -3 , ordinate is the ionization parameter, both 
in logarithm scale. Solid line is for log(CivA1549/AlinA1860), dot line is for 
log(SilvA1397/ Sim]A1892), dash line is for log(AllllA1860/ Sim] A1892), long dash 
line is for log(AllllA1860/SilvA1397), dash dot line is for log(ClvA1549/ Sim] A1892) 
and dash-triple-dot line is for log(SillA1814/ Sim] A1892). The point where the iso- 
contours cross determines the values of Logn e and LogU. The shaded area is the 
error bands for SillA1814. In 3C390.3 we add a long dash line for log(ClIl] A1909/ 
Sim] A1892). In J12014+0116 the shaded area is the error bands for SillA1814. See 
text for further details. 
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Fig. 12. — Correlation between rest frame equivalent width of AI111AI86O and 
Sim] A1892. Line in unweighted least square best fit. 
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Fig. 13. — Expected contribution from moderate density emitting gas as a function 
of density for ionization parameters log?7 = —2 (black lines) and \ogU = 2.5 (grey 
lines) . 
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Fig. 14. — Plane log U — log n in an expanded scale. The lower-left dot corresponds to 
a high density solution. Moderate density emitting gas is then added with increasing 
intensity of a Cm]A1909 component with the four Cm]A1909/ Sim] A1892 ratios 
reported in the figure. Connected positions corresponds to increasing density of 
Cin]A1909 emitting gas, logn = 9,9.5,10. Note that the case \ogn = 10 (grey 
symbols) is not appropriate for the observed data. 
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Fig. 15.— Contour plots for (a) J00103-0037, (b) J00521-1108, (c) J02390-0038, (d) 
J03036-0023, (e) J20497-0554 and (f) 23509-0052. Units and meaning of symbols 



are the same of Fig. 11 In the case of the objects shown in panels (c) and (f) the 
SiiiA1814 line is very weak and thus unreliable. For the objects in panels (a), (d) 
and (e) the line is also affected by the telluric absorption. For these objects we rely 
on the SilvA1397 line. 
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Fig. 16.— Contour plots for BAL quasars (a) J01225+1339 and (b) J01225+1339 
using Z0/A1304 an d (c) J01225+1339 using zcm]Ai909- Units and meaning of symbols 
are the same of Fig. [TTJ 
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Fig. 17. — Contour plots for corrected values considering the low density emission 
contribution of Clll]A1909. Contour plot for J00103-0037 is not shown because the 
correction is so large to be reliable. For the objects in panels (e) and (f) the SillA1814 
line is also affected by the telluric absorption. Abscissa, ordinate and symbols are 
the same as Fig. [TTJ 
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Fig. 18.— Contour plots for (a) J01225+1339, (b) J02287+0002 using z C ni]xi909, 
(c) J02390-0038, the extreme object (d) J12014+0116, and (e) J23509-0052, from 
the array of simulations computed for Z = 5Z & . Coordinates and symbols are as for 



Fig. 11 The intersection point improves in certain cases, but in others is the same. 
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Fig. 19. — Contour plots for the same objects of the previous figure from the array 
of simulations computed for Z = 5zT with ratios corrected because of low-density 
emission. Coordinates and symbols are as for Fig. 11 
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log Energy (Rydbergs) 

Fig. 20. — Spectral energy distribution used to compute the number of ionizing 
photons for Laor et al. (1997) in solid line and Mathews & Ferland (1987) in dotted 
line. Dashed line shows the Lyman limit. 
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Fig. 21. — M B h comparison for the high-z sample. Filled symbols refer to uncorrected 
intensity ratios; open symbols are for intensity ratios corrected because of low-density 
emission. Circles refer to solar metallicity; squares to 5 times solar metallicity and 
Si-Al enrichment. Abscissa and ordinate are logarithm of M B h in solar masses. Each 
point is labeled with the object name in short format. In ordinate we report Mbh 
values obtained with the method of this paper; in abscissa those obtained employing 
the Vestergaard & Peterson relationship described in the text. The shaded bands 
limit the 2a confidence level spread expected on the basis of the Vestergaard & 
Peterson relationship. 



